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Arithmetic  codes  are  useful  for  error-control  in  digital  computation  as  well  as  in 
data  transmission.  These  codes  are  especially  suitable  for  checking  or  correcting 
errors  in  arithmetic  processors  due  to  carry  propagation.  Two  known  classes  of 
arithmetic  codes  are  the  small-distance  high-rate  perfect  single-error  correcting 
codes  and  the  large-distance  low-rate  Man da lb aum-B arrows  codes.  These  codes  are 
analogous  to  the  Hamming  codes  and  the  maximum- length  sequence  codes  in  parity-check 
block  codes  respectively.  Most  other  arithmetic  codes  known  have  been  obtained 
by  computer-search.  The  discovery  for  a  systematic  way  of  constructing  arithmetic 
codes  with  intermediate-rate  and  intermediate-distance  has  been  the  sub-ect  of 
research  for  many  years.  Finding  simpler  decoding  algorithms  is  another  major 
unsolved  problem  in  arithmetic  codes.  Decoding  for  arithmetic  codes  by  matching 
the  orbits  or  permuting  the  residues  associated  with  the  codes  is  straightforward 
but  largely  impractical.  A  particularly  interesting  question  is  the  possibility  of 
decoding  arithmetic  codes  by  majority-logic.  In  this-thesis ,  we-h*ve  constructed  a 
class  of  intermediate- rate  intermediate-distance  binary  cyclic  arithmetic  codes.  A 
majority-logic  decoding  scheme  is  developed  for  the  code  constructed.  This  majority^ 
logic  decoding  scheme  is  also  applicable  to  a  large  class  of  cyclic  AN-codes  genex*atet 
by  the  primitive  cyclotomic  factors.  A  new  checking  technique  for  the  binary  adders 
has  been  developed.  This  separate  checking  technique  uses  less  redundancy  than  that 
required  for  b  triplicated  system,  and  the  decoding  procedure  is  simple.  The  appli¬ 
cation  of  majority-logic  decodab le  arithmetic  codes  to  the  error  control  in  high 
speed  multiplier  has  been  examined.  By  using  the  concepts  of  number  theory,  the  upper 
and  lower  bounds  -n  any  cyclic  AN~code  have  been  formulated. 
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Arithmetic  codes  are  useful  for  error-control  in  digital 
computation  as  well  as  in  data  transmission.  These  codes  are  especially 
suitable  for  checking  or  correcting  errors  in  arithmetic  processors  due 
to  carry  propagation. 

Two  known  classes  of  arithmetic  codes  are  the  small-distance 
high-rate  perfect  single-error  correcting  codes  and  the  large-distance 
low-rate  Mandelbaum- Barrows  codes.  These  codes  are  analogous  to  the 
Hamming  codes  and  the  maximum- length  sequence  codes  in  parity-check 
block  codes  respectively.  Most  other  arithmetic  codes  known  have  been 
obtained  by  computer- search.  The  discovery  for  a  systematic  way  of 
constructing  arithmetic  codes  with  intermediate-rate  and  intermediate- 
distance  has  been  the  subject  of  research  for  many  years. 

Finding  simpler  decoding  algorithms  is  another  major  unsolved 
problem  in  arithmetic  codes.  Decoding  for  arithmetic  codes  by  matching 
the  orbits  or  permuting  the  residues  associated  with  the  codes  is 
straightforward  but  largely  impractical.  A  particularly  interesting 
question  is  the  possibility  of  decoding  arithmetic  codes  by  majority- 
logic. 

In  this  thesis,  we  have  constructed  a  class  of  intermediate- 
rate  intermediate-distance  binary  cyclic  arithmetic  codes.  A  majority- 
logic  decoding  scheme  is  developed  for  the  code  constructed.  This 
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majority- logic  decoding  scheme  is  also  applicable  to  a  large  class  of 
cyclic  AN-codes  generated  by  the  primitive  cyclotomic  factors.  A  new 
checking  technique  for  the  binary  adders  has  been  developed.  This 
separate  checking  technique  uses  less  redundancy  than  that  required  for 
a  triplicated  system,  and  the  decoding  procedure  is  simple.  The 
application  of  majority- logic  decodable  arithmetic  codes  to  the  error 
control  in  high  speed  multiplier  has  been  examined.  Furthermore,  some 
theoretical  results  about  the  structure  of  arithmetic  codes  have  been 
obtained.  By  using  the  concepts  of  number  theory,  the  upper  and  lower 
bounds  on  any  cyclic  AN-code  have  been  formulated. 

In  summary,  the  results  presented  in  this  thesis  give  a  better 
belief  that  the  use  of  arithmetic  coding  for  the  error  control  in  digital 
computing  systems  is  indeed  promising. 
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1.  INTRODUCTION 


1.1  Error  Control  in  Digital  Computing  Systems 

The  rapid  growth  in  the  size  and  speed  of  modern  day  digital  comput¬ 
ing  systems  has  placed  stringent  reliability  demands  on  the  central  processors, 
especially  those  concerned  with  space  missions,  and  real  time  computer 
applications.  In  a  space  mission,  the  computers  require  high  reliability 
because  of  the  long  mission  time  during  which  manual  repair  of  a  failure  is 
impossible;  UK’s  a  real  time  computer  applications  require  high  reliability 
because  of  the  sori'us  consequences  error  might  cause  in  terms  of  inconvenience 
or  costly  mistakes.  Consequently,  the  error  control  techniques  that  can  be 
used  to  improve  the  reliability  of  a  digital  computing  system  are  very 
important. 

There  are  tx*o  major  techniques  to  deal  with  the  faults  in  a  digital 
system,  one  is  the  use  of  "software  diagnosis,"  and  the  other  is  the  design 
of  "hardware  checking  circuits."  In  the  software  fault  diagnosis  approach, 
a  specially  organized  checking  algorithm  is  incorporated  into  the  program 
and  the  arithmetic  results  are  verified  at  periodic  intervals  by  means  of 
additional  instructions  which  are  redundant  in  a  normally  functioning 
computer.  Correction  procedures  in  case  of  an  error  also  must  be  programmed. 
The  cost  of  this  approach  consists  of  the  additional  programming  effort,  the 
execution  time  and  storage  requirements  for  the  instructions.  The  logic 
design  of  the  arithmetic  processor  remains  unchanged.  An  alternative  to  the 
software  fault  diagnosis  is  the  design  of  hardware  checking  circuits.  In 
this  approach,  the  results  and/or  the  operands  of  each  arithmetic  operation 
arc  automatically  tested  for  acceptability  without  programmed  commands.  The 


indication  of  an  unacceptable  result  initiates  the  error  correction. 

The  cost  of  the  hardware  checking  approach  consists  of  the  additional 
logic  circuits  and  the  increased  execution  time  of  the  operations. 

Most  of  the  past  emphasis  to  improve  the  reliability  of  digital 
systems  has  been  on  making  more  careful  designs  and  using  more  reliable 
parts.  People  have  realized  that  it  is  also  possible  to  design  the 
system  that  it  operates  correctly  even  when  some  of  its  parts  fail. 

Such  approaches  invariably  use  redundancy.  The  above  mentioned  "soft¬ 
ware  diagnosis"  and  "hardware  checking"  techniques  are  both  the  approaches 
that  use  redundancy.  In  general,  by  utilizing  different  types  of  redun¬ 
dancy  or  inserting  the  redundancy  in  different  parts  of  a  system,  we  can. 
have  different  error  control  schemes  to  improve  the  reliability  of  digital 
computing  systems. 

The  most  commonly  used  scheme  for  hardware  checking  is  the  well- 
known  replication  of  processors.  Triplication  of  the  processor  with 
majority  output  voting  gives  error  correction;  and  duplication  of  the 
processor  with  output  comparison,  followed  by  further  diagnosis  in  the 
case  of  disagreement  gives  error  detection.  Both  methods  achieve 
checking  at  the  cost  of  complete  replication  of  the  processor.  An 
alternative  to  the  replication  of  processors  is  the  use  of  coding. 

The  encoding  of  the  operands  followed  by  the  application  of  decoding 
algorithm  to  the  results  requires  only  a  fractional  increase  in  cost, 
but  does  not  guarantee  complete  error  control. 

The  coding  approach  is  quite  different  from  the  technique  of  repli¬ 
cation  of  processors.  In  a  hardware  replication  scheme,  we  simply  organize 
a  number  of  less  reliable  components  to  form  a  reliable  computing  device; 
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while  in  a  coding  scheme  the  computing  device  is  fixed.  We  incorporate 
redundancy  into  data  themselves  which  are  being  processed;  in  other 
words  the  input  data  are  encoded  by  some  coding  technique  sc  that  the 
computational  errors  in  the  output  can  be  detected  and/or  corrected. 

For  example ,  if  we  use  an  AN-code,  the  operand  N  is  encoded  by  multi¬ 
plication  with  a  fixed  integer  A  before  being  represented  in  the  digital 
computer.  Thus,  only  a  fraction  1/A  of  the  possible  representations 
are  ever  intentionally  used  in  the  computer,  the  resultant  redundancy 
can  then  be  used  for  error  correction  and  detection. 

If  coding  technique  is  applied  to  improve  the  reliability  of  a 
computing  system,  there  is  no  need  to  supply  any  .software  diagnostic 
routines,  and  the  decoding  circuits  can  detect  and/or  correct  error 
immediately  after  it  occurs  so  that  any  unnecessary  propagation  of 
errors  can  be  avoided.  Also,  the  verification  of  the  acceptability  of 
input  operands  supplied  by  another  processor  provides  an  approach  to 
error-free  operation  of  systems  employing  intercommunicating  processors. 

Iii  fact,  some  of  the  computing  systems  designed  for  the  purpose  of 
being  used  in  spacecraft  on  long  missions  to  the  outer  planets  of  the 
solar  system  (e.g.  the  JPL  Self-Testing-and-Repair  Computer)  have 
employed  the  coding  scheme  to  fulfill  ultra-reliable  demands. 

Since  the  appearance  of  Shannon’s  theory  on  data  transmission,  the 
problem  of  how  to  design  an  efficient  scheme  by  which  the  information  can  be 
transmitted  reliably  across  noisy  channels  has  become  a  subject  of  continuing 
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importance.  The  invention  of  many  intelligent  codes  by  using  some  algebraic 
structures  has  highlighted  the  practical  aspects  of  the  coding  theory.  How¬ 
ever,  the  application  of  error-detecting  and  error-correcting  codes  to  a 
digital  computing  system  is  not  as  straightforward  as  it  has  been  to  a  data 
communication  system  such  as  space  coranunication  or  long-distance  telephone 
transmission.  Errors  in  a  digital  computing  system  are  of  different  nature 
from  the  errors  in  a  communication  channel.  Those  codes  which  perform 
powerfully  in  data  transmission  cannot  be  used  to  control  errors  caused 
oy  logic  faults  in  arithmetic  operations.  In  general,  there  are  two  types 
of,  errors  that  occur  in  a  digital  computing  system,  transmission  errors 
and  arithmetic  errors.  Transmission  errors  usually  occur  in  a  large 
digital  system  when  many  intercommunicating  data  processors  are  employed, 
the  data  flow  across  different  units  may  be  affected  by  noise.  These  errors 
are  of  the  same  nature  as  errors  in  a  communication  system.  Arithmetic 
errors  occur  in  the  operation  of  an  arithmetic  processor  whenever  its 
actual  behavior  deviates  from  the  expected  behavior.  Those  errors  are 
caused  by  faults  of  logic  circuits  within  the  processor.  The  main  difference 
between  transmission  errors  and  arithmetic  errors  is  that  the  latter  “ay 
cause  carry  or  borrow  propagations  so  that  a  single  fault  in  logic  may 
cause  several  output  bits  to  be  in  error.  Therefore,  an  efficient  error¬ 
detecting  or  error-correcting  code  used  in  a  digital  computing  system 
must  be  able  to  deal  with  both  types  of  errors.  The  codes  that  are 
specially  designed  for  the  use  of  detecting  and/or  correcting  errors  in 
digital  computation  as  well  as  data  transmission  are  called  arithmetic 
codes.  It  is  the  purpose  of  this  thesis  to  study  the  error  control  in 
digital  computing  systems  by  using  arithmetic  coding. 
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The  theory  of  arithmetic  coding  has  drawn  the  attention  of  many 
researchers  since  Diamond  [20]  first  proposed  his  paper  in  1955.  Much 
research  has  been  done  on  detecting  and  correcting  errors  in  arithmetic 
operations.  Most  of  this  work  has  been  concerned  with  the  arithmetic  AN-code, 
(we  shall  make  somewhat  detailed  discussion  on  the  fundamental  concepts  and 
definitions  of  arithmetic  AN-codes  in  the  following  chapters.)  Massey  [36] 
presented  an  excellent  survey  of  the  early  work  on  these  codes.  The 
popularity  of  these  codes  is  due  to  the  fact  that  they  have  a  nice  algebraic 
structure  and  can  be  used  in  a  computer  with  little  or  no  change  in  the  cir¬ 
cuit  design.  Although  many  important  results  have  been  obtained  in  the  theory 
of  arithmetic  codes,  there  are  still  two  major  unsolved  problems: 

(1)  The  problem  of  how  to  find  a  systematic  way  to  construct  a  nultiple 
error-correcting  arithmetic  code.  The  existing  codes  are  either  of 
high-rate,  low  error-correcting  capability  or  of  high  error  control 
power  but  low  information  rate.  The  codes  with  intermediate  rate 
and  error  control  capability  are  usually  obtained  by  computer  search. 
Therefore,  a  systematic  way  of  constructing  good  arithmetic  AN-code 
is  much  in  need. 

(2)  A  more  important  problem  is  the  implementation  (or  decoding)  of  the 
arithmetic  codes.  The  practicality  of  using  coding  in  a  digital 
system  depends  heavily  on  how  complex  the  decoders  are.  Unfortunately, 
the  decoding  problem  of  arithmetic  codes  is  quite  complex  in  general 
because  of  the  carry  propagation  caused  by  errors.  The  only  known 
decoding  method  is  the  "permutation  of  residues"  [29,30],  which  is 
basically  a  table-look-up  approach  and  impractical.  How  to  find  a 


decoding  scheme  that  can  be  easily  implemented  is  then  a  problem  of 
particular  significance.  One  approach  is  suggested  by  the  following 
consideration:  since  the  majority-logic-decoding  (MLD)  scheme 
[6]  is  by  far  the  easiest  one  that  can  be  applied  in  a  variety  of 
ways  to  communication  channels,  it  would  certainly  be  of  interest  to 
see  if  the  decoding  of  arithmetic  codes  could  also  be  handled  by 
using  majority  logic. 

The  solutions  of  these  two  problems  will  greatly  increase  the 
feasibility  and  effectiveness  of  using  arithmetic  coding  scheme  in  digital 
computing  systems.  These  are  two  of  the  topics  that  are  investigated  in 
this  thesis. 

1.2  Scope  of  the  Thesis 

The  main  objectives  of  this  thesis  are  to  construct  new  arithmetic 
codes  with  better  performance,  to  devise  decoding  schemes  for  arithmetic 
codes  that  can  be  practically  implemented,  to  develop  efficient  checking 
techniques  and  checking  model  for  arithmetic  operations  in  computer,  and 
to  study  the  theoretical  structures  of  the  arithmetic  codes. 

Chapter  2  is  a  general  review  of  the  background,  definitions  and 
concepts  that  are  helpful  to  the  discussion  of  the  material  presented  in 
later  chapters.  Concepts  of  error  weight  and  code  distance,  used  for 
determining  the  codes'  error  correcting  capability  are  discussed.  The 
cyclicity  of  the  codes  and  the  analogies  between  the  AN-codes  and  parity- 
check  codes  are  briefly  examined. 

Chapter  3  considers  tho  construction  of  new  arithmetic  codes  with 
multiple-error-correcting  capability.  A  fast  algorithm  for  finding  the 
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binary  representation  forms  of  certain  integers  is  derived.  By  using  the 
number  theoretical  properties  of  these  integers,  a  class  of  new  codes  is 
constructed.  The  error  control  power  of  these  codes  be  determined 
analytically.  The  root-distance  relation  of  a  class  of  cyclic  AN-code 
is  investigated,  the  result  is  completely  analogous  to  the  BCH  bound  for 
the  cyclic  parity-check  codes.  In  addition,  the  bounds  on  the  error- 
correcting  capability  of  any  cyclic  AN-codes  are  determined. 

Chapter  4  describes  a  majority- logic  decoding  scheme  for  some  cyclic 
AN-codes.  This  decoding  scheme  is  completely  different  from  the  permutation 
of  residue  approach,  the  implementation  of  the  codes  can  be  achieved  by 
using  majority  logic  alone .  The  application  of  the  majority-logic  decoding 
algorithm  to  the  cyclic  AN-codes  constructed  in  Chapter  3  is  examined,  and 
a  multi-step  majority  decoding  scheme  that  can  be  applied  to  larger  classes 
of  cyclic  AN-codes  is  obtained. 

Chapter  5  suggests  a  new  spearate  concurrent  checking  technique  for 
the  arithmetic  operations  in  a  digital  computer.  The  application  of  error- 
correcting  AN-codes  shows  the  effectiveness  of  the  separate  checking  tech¬ 
nique.  Finally,  a  study  jf  the  error  control  in  high  speed  arithmetic  is 
included. 

Some  conclusions  on  the  results  presented  in  this  thesis  and  sugges¬ 
tions  on  possible  future  work  are  also  discussed  in  the  last  chapter. 
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2.  PRELIMINARIES  AND  BACKGROUND 

2.1  The  Arithmetic  AN-Code 

The  arithmetic  AN-code  is  specially  designed  for  the  use  of  correcting 
and/or  detecting  errors  in  digital  computation  as  well  as  data  transmission. 

The  code  is  a  set  of  integers  of  the  form  AN,  where  A  is  a  fixed  integer  called 
the  generator  of  the  code.  Fcr  each  integer  N  from  the  information  set 
Zj  =  {0,1,2,.,.,  B-l],  there  is  a  corresponding  code  number  AN,  hence  the  integer 
N  is  "encoded"  by  multiplication  with  the  fixed  generator  A  before  being  repre¬ 
sented  in  the  computer  arithmetic.  Thus  only  a  fraction  1/A  of  the  possible 
representations  are  actually  used  in  the  computer,  the  resultant  redundancy 
can  then  be  used  for  error  detection  and  correction. 

The  AN-code  possesses  linear  property,  for 

ANj  +  AN2  *  A(Nj+N2) 

so  that  the  coded  form  for  the  sum  of  two  information  integers  is  the  sum  of 
the  coded  integers.  Therefore,  two  coded  numbers  can  be  added  in  an  ordinary 
adder;  if  the  sum  is  not  a  properly  coded  number  error  is  checked. 

As  mentioned  in  the  introductory  chapter,  the  errors  that  occur  in  a 
digital  computing  system  are  more  complicated  than  the  transmission  type  of 
errors.  Consequently,  the  identification  of  the  likely  modes  of  failore  in  a 
computing  system  and  the  formulation  of  a  convenient  description  of  the  resultant 
errors  in  its  output  must  be  considered  differently  from  that  of  a  data  trans¬ 
mission  system. 
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The  measure  of  the  weight  of  an  arithmetic  error  and  the  error- 
correcting  capability  of  an  AN-code  have  been  developed  by  several  authors 
[35,37],  a  brief  summary  of  these  results  will  be  made  in  the  following  sections. 


The  circuitry  for  adding  two  integers  in  radix-2  form  can  be  con¬ 
sidered  as  a  sequence  of  elemental  adding  units.  Suppose  that  two  binary  inte¬ 
gers  I-  *  (a  a  ,...a,an)  and  I  »  (b  b  ,...b,b_)  are  added,  then  the  i-th 
1  n  n-1  1  O'  2  n  n-1  1  0' 

elemental  adding  unit  performs  the  addition  of  ai+^,  b^+^  .and  the  carry  bit 
from  the  (i-l)-th  unit.  Each  unit  orms  a  sum  bit  and  a  carry  bit  to  the  next 
unit.  A  "single  failure"  in  the  adder  is  then  considered  as  an  incorrect  sum 
formed  by  one  adding  unit  or  the  generation  of  an  incorrect-carry  bit  by  one 
adding  unit;  hence  a  single  failure  the  i-th  unit  may  cause  an  error  of 
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value  +  2*  (due  to  incorrect  sum)  or  +  (due  to  incorrect  carry) .  If  the 

erroneous  sum  differs  from  the  correct  sum  by  an  amount  E,  it  is  then  quite 
natural  to  think  of  the  weight  of  E,  i.e.  the  number  of  errors,  as  the  least 
number  of  terms  of  the  form  +  2*  or  -2^  whose  sum  is  E.  This  measure  of  the 
weight  of  an  error  is  called  "arithmetic  weight"  and  can  be  formalized  as 
follows: 

Definition  2.2.1:  The  arithmetic  weight  of  an  integer  E  (which  may  be 
positive,  negative  or  zero),  denoted  W(E),  is  the  minimum  number  of  non-zero 
coefficients  in  the  modified  binary  form  of  E: 

E  «  e^  +  e^  2^  +  2^  +  ...+  e^  2^  +  ... 

where  (2.2.1) 

■  0,  1,  or  -1 

Although  the  modified  binary  form  of  an  integer  is  not-  unique  in 
general,  there  is  a  particular  modified  binary  form  introduced  by  Reitwiesner 
[43]  which  is  unique  and  which  has  been  proved  of  much  importance  in  the  theory 
of  arithmetic  code,  called  the  non-adjacent-form  (or  NAF  for  short).  A  modi¬ 
fied  binary  form  is  said  to  be  a  NAF  if  the  coefficients  e^'s  in  (2.2.1)  satisfy 
ei  ei+l  "  0  for  i  =  0,  1,  2,  ...  ,  in  other  words  if  there  are  no  two  adjacent 
non-zero  coefficients  in  the  form. 

The  most  important  property  of  the  NAF  of  an  integer  I  is  that  no 
modified  binary  form  of  I  has  fewer  non-zero  coefficients  than  the.  NAF  of  I  [43] . 
Therefore,  the  arithmetic  weight  of  any  integer  equals  the  number  of  non-zero 
coefficients  in  its  NAF. 


11 


Definition  2.2.2:  The  arithmetic  distance  between  the  integers  1^  and 
denoted  is  the  arithmetic  weight  of  their  difference,  i.e.  D^,^) 

-  wo^-ip. 

It  can  be  shown  [35]  that  the  arithmetic  distance  satisfies  the 
following  three  properties: 

DO^.ip  -  D(I2,I1)  (symmetry) 

D(Ij^I2)  Ss  0  with  equaltiy  iff  1^  =  I2  (positive  definite) 
D(I^,I2>  S  D(IltI3)  +  D(I3,I2)  (triangle  inequality) 

Hence,  any  set  of  integers  with  arithmetic  distance  taken  as  the  measure  of 
"distance"  fora  a  metric  space. 

Definition  2.2.3:  The  minimum  arithmetic  distance,  D  ,  ,  of  an  arithmetic 

min 

AN-code  is  the  minimum  of  the  arithmetic  distances  between  all  pairs  of  distinct 
code  words  in  the  code. 

Since  the  difference  of  two  code  words  is  another  code  word,  the  D  , 

’  min 

of  an  arithmetic  AN-code  is  equal  to  the  minimum  arithmetic  weight,  ^ >  of  the 
B-l  non-zero  code  words  in  the  code. 

The  error-correcting  capability  of  an  AN-code  is  completely  deter¬ 
mined  by  the  minimum  distance  of  the  code  as  Massey  [35]  has  shown  that  an 

AN-code  with  D  .  =  2t  +  1  can  correct  all  arithmetic  errors  of  weight  t  or 

min 

less,  and  an  AN-code  with  D  .  =  t  +  1  can  detect  all  arithmetic  errors  of 

min 

weight  t  or  less.  Therefore,  the  knowledge  of  minimum  distance  of  a  given 
arithmetic  AN-code  has  become  a  very  important  subject. 
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2.3  Cyclic  AN-Code  and  Modular  Distance 


Suppose  that  an  integer  I  has  a  binary,  n-tuple  form  (a^  1  a^^  . . . 


Sq) .  We  denote  by  P(I)  the  integer  whose  binary  n-tuple  form  is  the  left  cyclic 

shift  of  that  for  I,  namely  (a  .a  ,  a.  a  .).  If  the  largest  code  word 

J  n-2  n-3  0  n-1 

A(B-l)  in  an  AN-code  requires  n  bits  for  its  radix-2  form,  then  we  define: 


Definition  2.3.1:  An  AN-code  is  cyclic  if  its  set  of  codewords  is  closed 
under  cyclic  shifting,  i.e.  if  for  every  code  word  AH,  the  integer  P(AN)  is 
another  code  word. 

It  has  been  shown  [11,36]  that  the  generator  A  and  the  number  of 

codewords  B  in  a  cyclic  AN-code  satisfy  AB  *  2e^^  -  1,  where  e(A)  denotes  the 

least  positive  integer  such  that  A  divides  26^  -1,  and  is  called  the  exponent 

of  2  modulo  A  [49].  Therefore,  the  block  length  of  a  cyclic  AN-code  generated 

by  A  is  e(A).  Mathematically,  ^’l  the  codewords  in  a  cyclic  AN-code  form  a 

efAl 

principal  ideal  in  the  ring  of  integers  modulo  2  -1,  with  A  as  the  ideal 

generator. 

0  (A) 

In  a  cyclic  AN-code,  the  addition  of  codewords  is  the  modulo  2  -1 

arithmetic,  a  negative  integer  -I  is  represented  by  its  one's  complement,  i.e. 

0  (A) 

2  -1-1.  The  arithmetic  weight  of  negative  integers  then  is  defined  for 

their  one's  complements.  This  leads  naturally  to  a  modified  measure  of  the 
weight  of  the  integers  in  the  ring  of  integers  modulo  2n-l: 


Definition  2.3.1:  The  "modular  weight"  of  an  integer  F  in  the  ring  of 
integers  modulo  2n-l  is  the  minimum  of  W(F)  and  W(2n-1-F),  and  is  denoted 

Wm(F> • 

m 


-zswTrzrju 
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Definition  2.3.2:  The  "modular  distance"  between  the  integers  1^  and 

in  the  ring  of  integers  modulo  2n-l  is  the  modular  weight  of  their  difference 

and  is  denoted  D  (I.,I,),  i.e.  W  (1,-1.). 

mi/  ml/ 

With  the  definitions  of  modular  weight  and  modular  distance,  we  can 
show  that  the  error  correcting  capability  of  a  cyclic  AN-code  in  modulo  arith¬ 
metic  is  completely  determined  by  the  modular  minimum  distance  of  the  code,  a 
result  which  is  analogous  to  the  fundamental  result  stated  in  the  previous 
section. 

2.4  Arithmetic  Codes  for  Single  Error 

Single  error  detecting  and  single  error  correcting  arithmetic  AN-codes 
have  been  thoroughly  investigated  by  Brown  [8]  and  Peterson  [37].  Some  well- 
krown  results  will  be  summarized  in  the  following  theorems. 

Theorem  2.4.1:  An  AN-code  with  generator  A  >  1  and  odd  can  detect  any 
single  ertor  in  radix-2  arithmetic. 

Theorem  2  4.2:  Let  M2(A,d)  be  the  smallesc  positive  integer  whose  product 
with  A  has  arithmetic  weight  less  than  d,  then  the  arithmetic  AN-code  with 
0  <  N  <  M^A.d),  has  Dmin  4  d. 

Theorem  2.4.3:  Let  A  >  1  be  an  odd  integer,  e  be  the  smallest  integer  that 
satisfies  2e  s  1  mod  A,  then  >^(A,3)  *  (2e-l)/A.  Similarly,  if  2e 
then  M2  (A, 3)  =  (2e+l)/A. 


s  -1  mod  A, 


14 


Theorem  2.4.4:  If  A  is  an  odd  prime  and  if  2  is  a  primitive  root  of  A, 
e 

then  M^(A;3)  =  (2^  +1)/A.  Similarly,  if  -2  is  the  primitive  root  of  A,  then 
e 

^2^,3)  =  (22  -1)/A,  where  e  =  A-l. 

The  above  results  yield  a  way  of  constructing  any  single  error-detect¬ 
ing  and  single  error-correcting  arithmetic  AN-codes.  More  importantly,  the  single 
error-correcting  codes  mentioned  above  are  perfect  codes  or  sphere-packed  codes, 
i'-e.,  any  integer  in  the  ring  of  integers  modulo  2e  -1  is  at  most  distance-one 
away  from  some  code  word.  These  codes  are  analogous  to  the  well-known  Hamming 
codes  in  parity-check  block  codes. 

2.5  Mandelbaum-Barrows  Equidistant  Codes 

The  Mandelbaum-Barrows  codes  [5,32]  are  the  first  systematically  con¬ 
structed  class  of  AN-codes  with  >  3.  They  are  also  the  first  AN-codes  to 
be  recognized  as  cyclic.  The  structure  of  these  codes  is  simple: 

Theorem  2.5.1:  Let  B  be  :an  odd  prime  with  2  as  its  primitive  root,  let  the 
generator  A  of  the  cyclic  AN-code  be: 

A  =  (2e(B)-l)/B  =  (2B_1-1)/B 

Then  the  minimum  distance  of  the  code  is  T (B4-l)/31  ,  where  Ixl  denotes  the  inte¬ 
gral  part  of  x> 

For  a  properly  chosen  B,  these  codes  can  correct  any  multiple  arith¬ 
metic  errors.  Moreover,  these  codes  also  possess  an  interesting  property  that 
all  the  nonzero  code  words  are  some  cyclic  shifts  of  the  generator  A,  hence  they 
all  have  the  same  arithmetic  weight  f(B+l)/31  .  These  codes  are  analogous  to  the 
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maximum- length  sequence  codes  in  parity-check  block  codes,  they  are  of  large- 
distance  but  low-rate,  and  stand  at  the  opposite  end  of  the  coding  spectrum  from 
the  perfect  single-error-correcting  codes  of  the  previous  section. 

2.6  Analogy  Between  Arithmetic  AN-Codes  and  Parity-Check  Codes 

We  begin  with  a  brief  review  of  the  theory  of  cyclic  parity-check  codes 
Extensive  treatments  of  this  subject  can  be  found  in  Peterson  [37]  and  Berlekamp 
[6]. 

n-i  t 

With  a  polynominal  f(x)  =2  f  x  of  degree  less  than  n  with  coefficients 

i=o  1  ’ 

in  a  field  F,  we  associate  the  vector  f=[f  .,  f  _,...,  f . ,  fn]  in  the  vector 

n-1  n-2  1  0J 

space  Fn.  A  parity-check  code  is  simply  a  set  of  such  vectors  which  form  a  sub¬ 
space  of  the  linear  vector  space  Fn.  The  Hamming  weight  of  a  vector  (or  a  code 
word)  is  defined  as  the  number  of  non-zero  components  in  the  vector,  and  the 
Hamming  distance  of  two  vectors  is  the  number  of  positions  in  which  the  two 
vectors  differ.  Since  the  vectors  are  closed  under  addition,  the  minimum  Hamming 
distance,  d  ,  °f  such  a  code  equals  the  minimum  Hamming  weight  of  the  non-zero 
vectors  in  the  code. 

A  parity-check  code  is  cyclic  if  the  cyclic  shift  of  every  codeword  is 
also  a  code  word.  If  we  denote  the  cyclic  shift  of  f  by  f ' ,  then  the  correspond¬ 
ing  polynomial  f '  (x)  of  V  is  obtained  by  multiplying  f(x)  by  x  and  take  the 
residue  modulo  xn-l.  It  can  be  shown  [37]  that  the  code  words  of  a  cyclic  parity 
check  code  form  a  principal  ideal  in  the  polynomial  ring  modulo  xn-l  over  some 
field  F.  Therefore,  each  code  word  has  a  corresponding  polynomial  which  is  a 
multiple  of  some  fixed  polynomial  g(x),  called  the  generator  polynomial  of  the 
ideal.  In  addition,  g(x)  is  monic  and  divides  xn-l.  The  polynomial 
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h(x)  =  (xn-l)/g(x)  is  called  the  parity-check  polynomial  of  the  code,  and  n,~ 
the  length  of  the  code  is  chosen  as  the  least  positive  integer  such  that  g(x) 
divides  xn-l. 

In  light  of  the  above,  the  analogy  of  cyclic  parity-check  code.,-  ith  the 
cyclic  arithmetic  AN-codes  it.  evident.  The  analogous  quantities  are: 


Cyclic  parity-check  codes 
x 


g(x) 

h(x) 


period  of  g(x) 

d  . 
min 


Cyclic  arithmetic  AN-codes 
2 
A 
B 

2n-l 

exponent  of  2  mod  A 

D  . 
min 


The  closeness  of  these  analogies  between  cyclic  parity-check  codes  and 
cyclic  AK-codes  strongly  implies  that  there  exist  cyclic  AN-codes  analogous  to 
many  important  cyclic  parity-check  codes,  in  particular,  the  Bose-Chandhuri- 
Hocquenghem  codes  and  the  majority- logic  decodable  geometry  codes.  Many  investi¬ 
gations  have  been  made  on  this  problem,  but  to  date  very  little  is  known. 

2 . 7  Remarks 

There  are  other  important  results  and  concepts  in  the  theory  of  arith¬ 
metic  codes  that  have  not  been  mentioned.  It  should  be  remarked  that  the  pur¬ 
pose  of  this  chapter  is  to  provide  some  of  the  terminology  and  concepts  which  will 
be  used  with  regard  to  the  following  discussion.  The  material  presented  here  is 
just  a  brief  review  of  the  arithmetic  coding  theory,  more  detailed  treatments  of 
the  related  subjects  can  be  found  in  [35,36,37]  etc. 
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3.  CONSTRUCTION  OF  MULTIPLE-ERROR-CORRECTING 


ARITHMETIC  AN-CODES 


3.1.  Introduction 

The  theory  of  arithmetic  coding  has  been  developed  since 
L'^moad  [20]  first  proposed  his  paper  in  1555.  Single  error-correcting 
arithmetic  codes  have  been  investigated  by  Brown  [8  3,  Peterson  [83], 
and  Bernstein  [7  ].  Some  results  of  multiple  error-correcting  arithmetic 
codes  have  been  reported  by  Barrows  [ 5  ] »  Mandelbcum  [ 32] ,  Chang  and 
Tsao-Wu  [ll],  Chien,  Hong  and  Preparata  [  18] . 

Since  the  error-ccrrecting  capability  of  arithmetic  codes  is 
directly  related  to  the  minimum  distance  of  ^he  codes,  an  analytical  way 
of  calculating  the  distance  for  arithmetic  codes  is  important.  Chien, 
Hong  and  Preparata  [13],  Tsao-Wu  and  Chang  [  11]  independently  discovered 
a  computational  algorithm  for  minimum  distance  of  cyclic  AN-codes.  By 
using  number  theoretic  concepts,  they  divide  the  code  words  of  a  cyclic 
AN-code  into  a  number  of  disjoint  sets,  each  set  of  code  words  is  called 
an  orbit.  The  code  words  in  the  same  orbit  are  the  AN's  such  that 
N  =  (k«2^)  mod  B,  for  some  fixed  k,  and  i  =  0,1,2,. ..,n,  AB  =  2n-l. 
Therefore,  each  k  defines  an  orbit.  It  can  be  seen  that  the  code  words 
in  the  same  orbit  are  of  equal  weights,  hence  the  minimum  weight  of  the 
code  is  the  minimum  of  all  the  weights  of  different  orbits. 

Perhaps  the  most  interesting  property  of  this  orbital  theory 
is  the  following  theorem: 

Theorem:  The  weight  of  a  code  word  AN,  represented  in  binary  form 


AN  =  a  ,a  «,.... a, a» 
n-1  n-2  1  0 
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where  afl_^  ■  ((N«2*)mod  B)  mod  2,  1  <  i  <  n,  is  given  by  the  number  of 
residues  N*2^  mod  B€M^;  where  M^,  the  middle-third  region  of  B  is  defined 

Mfi  =  [x|b  <  3x  <  2B} . 

The  result  of  this  theorem  enables  us  to  calculate  the  arithmetic 
weight  of  each  orbit  by  counting  the  number  of  residues  in  M„.  Computation" 

D 

ally,  this  orbital  algorithm  might  still  be  difficult  especially  when  the 
code  contains  a  large  number  of  orbits.  An  easier  algorithm  to  find  the 
minimum  distance,  or  perhaps  some  closed-form  formula,  is  necessary.  Further¬ 
more,  in  practical  use  the  discovery  of  systematic  way  of  synthesizing  a 
code  which  corrects  a  specified  number  of  errors  is  more  important.  The 
problem  that  how  to  construct  a  code,  rather  than  analyze  a  given  code, 
is  therefore  of  both  theoretical  and  practical  significance. 

In  this  chapter,  we  shall  construct  a  class  of  cyclic  AN-codes 
whose  generators  are  the  products  of  cyclotomic  factors.  The  multiple 
error-correcting  capability  of  these  codes  can  be  determined  exactly. 

These  codes  fill  in  the  gap  between  the  single-error-correcting  codes 
and  the  Mandelbaum-Barrows  codes  as  far  as  the  distance  and  rate  are 
concerned.  The  construction  of  these  codes  is  based  on  some  interesting 
properties  possessed  by  the  binary  forms  of  certain  integers  which  are 
used  as  the  generators  of  the  codes. 

We  shaJ 1  first  establish  a  fast  algorithm  for  finding  the  binary 
forms  of  certain  integers,  by  using  the  complementary,  symmetric  properties 
we  shall  construct  a  class  of  cyclic  AN-codes.  Also,  the  analogy  between 
cyclic  AN-codes  and  cyclic  polynomial  codes  will  be  studied,  we  shall 


demonstrate  a  root-distance  relationship  of  a  class  of  AN-codes  which  is 
shown  to  be  exactly  the  same  as  the  well-known  BCH  bound  on  the  distance 
of  cyclic  polynomial  codes.  Furthermore,  by  using  number  theoretic 
concepts,  the  strict  upper  and  lower  bounds  on  the  error-control  capability 
of  any  cyclic  AN-codes  are  derived.  In  general,  the  results  obtained  in 
this  chapter  are  all  theoretical,  the  possible  practical  applications  and 
implementations  of  the  codes  constructed  will  be  discussed  in  the  following 
chapters. 


The  representation  of  integers  in  binary  form  has  been  under 

investigation  for  various  applications.  In  particular,  the  binary 

representation  is  of  importance  in  fast  computer  arithmetics,  number 

systems,  algebraic  coding  theory  and  arithmetic  codes.  In  this  section, 

nl 

the  binary  representation  of  the  integer  A  of  the  form  (2n-l)/(2  -1). 

n2  , 

(2  -1),  with  n^’n2  **  n  n^>n2  being  two  relatively  prime  integers, 

is  considered.  Firstly,  some  interesting  properties  of  the  binary 
representation  form  will  be  presented.  It  will  be  shown  that  the  binary 
form  of  A  has  l's  and  0's  at  certain  fixed  positions  with  complementary 
symmetrical  properties.  Secondly,  a  simple  and  fast  algorithm  for  finding 
the  binary  form  will  be  introduced.  From  this  algorithm,  the  binary  form 
of  a  large  integer  can  be  easily  obtained  without  any  calculation. 
Formulations  and  Notations 

Let  A  be  an  integer  of  the  form 

n  n2 

A  =  (2-1)/ (2  -1) (2  -1) 
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where  and  ng  are  two  relatively  prime  positive  integers  with  >  n^ 
and  n  «  n^n^.  The  integer  A  can  be  uniquely  expressed  in  the  following 
binary  form 

n-1 

A  *  [  a.21  with  a,  =  0  or  1  for  0  <  i  <  n-1 
i=0 

ni 

Let  B  ■  (2  -1)(2  -1),  then  AB  ■  2n-l.  The  coefficients  a^  in  tlw  binary 

representation  form  of  A  can  be  determined  [12]  by 


ii"  i 

a^  «  (2  mod  B)  mod  2 


(3.2.1) 


ii"i  ti"i 

where  (2  1  mod  B)  is  the  residue  of  2  modulo  B.  The  above  equation 

II*  i 

states  that  a  =  0  if  and  only  if  (2  mod  B)  is  an  even  integer. 

nl 

Consider  the  binary  form  of  the  integer  A* 2  , 

nl 

C  -  A  ♦  2 
n-1 


E  c  A 

i=0 


c^  -  0  or  1  for  0  <  i  <  n-1 


By  (3.2.1), 


,  n. 

c^  =  [2°  *2  mod  B]  mod  2 


(3.2.2) 


Similarly,  let 


nl 

D  -  A  •  (2  -1) 
n-1 

*  E  d  21 
i»0 


d^  =  0  or  1  for  0  <  i  <  n-1 


then, 
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Let  us  denote 

[2°  *  •  2  ^  mod  B]  *  Xj, 

r-.n-i  .  _ 

L2  mod  BJ  ■>  y± 

‘I 

[2n“1(2  -1)  mod  B]  -  z± 

where  0  <  XjL,  y ±,  zj[  <  B.  Note  that  c±  s  Xj[  mod  2,  =  y±  mod  2  and 

di  S  zi  mod  2* 

Properties 

Property  1:  C  ■  0  and  a  «  1. 

o  o 

Proof:  Since  B  divides  2n-l,  2nsl  mod  B  and  y  »  1.  In  addition, 
n  nl  nl  nl 

since  2  *2  s2  mod  B,  x  «  2  .  Thus  c  «  0  and  a  «  1. 

o  o  o 

Q.E.D. 

Suppose  that  i  >  0.  Let  n-i  «  k.  By  Euclidean  division 
algorithm  we  have  k  ■  qn^+r,  where  0  <  r  <  n^  Notice  that 
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nl  n2 

B  -  (2  -1) (2  -1)  implies: 


n«  ni 

2  Z(2  -1)  b  2  -1  mod  B 


Thus 


jn„  n.  (j-l)n„  n 

2  (2  -1)  s  2  (2  -1) 


n. 


=2  -1  mod  B  for  j  >  0 


Now,  equation  (3.2.5)  can  be  used  successively  to  obtain 


(3.2.5) 


.  .  k-n.  n,  k-n. 

2n-i  -  2k  o  2  X(2  *-l)  +2  1 

(k-n..)  mod  n  n.  k-2n..  n.  k-2n.. 

h  2  2(2  -1)  +  2  i(2  -1)  +  2 

n  (k-n.)  mod  n„  (k-2n.)  mod  n  k-3n  n1  k-3n. 
h  (2  -1)[2  2  +  2  Z]  +  2  L(2  -l)+2 

n..  (k-n.)  mod  n_  (k-2n.)  mod  n  (k-3n..)  mod  n. 

=  (2  -1)[2  1  2  +  2  1  2  +  2  1  2]  + 


k-4n1  n1  k-4n 

2  A(2  -1)  +  2  1 


n.  (k-n1 )  mod  n„  (k-2n. )  mod  n_  (k-qn1 )  mod  n0 

s  (2  -1)[2  2  +  2  2  +  ...  +  2  Z]  + 

k-qn. 

2  (mod  B) 


Let 


(k-n^)  mod  n ^ 


t  •  f 


S  »  [2 


(k-2n.)  mod  n„ 

+  2  +  •••  +  2 
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(k-qn^)  mod  02 


] 


Then 


n. 


2n"i  «  2k  e  (2  1-1)  S  +  2r  mod  B 


(3.2.6) 


Similarly, 


n,  n. 


,  ...  ...  n.  k  mod  n. 

2n-1  •  2  =  (2  -1)  Sf2r  +  (2  -1)2 


mod  B  (3.2.7) 


Notice  that  there  are  q  terms  in  S,  each  term  is  a  power  of  2.  In 
addition,  these  terms  are  distinct,  for  otherwise, 


k-mn^  =  k-m'n^  mod  n«  for  some  0  <  ra,  m'  <  q  < 

This  would  imply  m-m'sO  mod  n2.  Since  (n^,^)-!  and  0  <  m,mf<  n2»  the 

congruence  is  impossible  unless  m  ■  m'. 

For  i  >  0,  q  <  n  -1.  Since  all  the  q  terms  in  S  are  distinct, 
n,  1  «  n2"^  « 

q  <  ng,  and  2  -1  «  1+2  +2  +...+2  ,  we  have  S  <  2  -2.  Also, 

n1 

2r  <  2  -1;  thus: 


n  n.  n  n 

(2  -1)  Sf2r  <  (2  -1)  (2  -2)  +  (2  -1)  -  B 

-  j  o. 

By  (3.2.6),  the  residue  of  2  modulo  B  is  equal  to  (2  -1)  Sf2  . 

Thus  we  have 


Property  2: 


4  1 

y^^  «•  (2n*i  mod  B)  ■  (2  -1)  Sf2r,  where  i  ^  0 

(k-n1 )  mod  n  (k-2n. )  mod  n„  (k-qn. )  mod  n, 

S-2  1  +2  1  +...+  2 


■  n-i  ■ 


0  <  r  <  nx 


k 


qnL  +  r. 


(3.2.8) 


’v  i»»  ^ -jt n '.■r'«J  ■n'M*v  •jw?»r"»r 
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Similarly,  the  following  property  can  be  obtained: 

n2  ^  . 

Property  3:  y±  =  (2  -1)  T  +  2  ,  where  i  f  0 

(k-nj  mod  n.  (k-2n„)  mod  n.  (k-pO  mod  n 

T  *  2  2  +  2  2  1  +...+2 

k  »  n-i  =  p  *  n^  +  ti,  0  <  t  <  n2  (3.2.9) 


Theorem  3.2.1:  <  yt  for  0  <  i  <  n^ 

Proof:  For  0  <  i  <  n^  n  -  i  =  n^  -  i  *  (n2~l)n1+(n1-i) ,  we  have 


q  »  n2-l,  r  ■=  n^i 

_  j  n.  nino+ni “i  r  r 

Since  2  1  ♦  2  1  =  2  1  2  1  =212*2=2  mod  B, 


4  nl 

x.  *  (2  •  2  mod  B)  -  2 


(3.2.10) 


By  (3.2.8),  we  have  y±  >  x±. 


Q.E.D. 


Property  4:  c±  »  0,  ai  =  1  fc-  0  <  i  <  n^  and  c  -1,  an^=  0. 

Proof:  By  (3.2.10),  x±  is  even  if  0  <  i  <  ^  and  xn  is  odd. 

Therefore,  c.  **  0  for  0  <  i  <  n.  and  c  =1.  From  theorem  3.2.1, 
i 

xt  -  y±  <  0.  Note  that  z±  =  x1~yi  mod  B,  we  have 

zi  “  B  +  Xi  ’  yi 


From  (3.2.4),  we  see  that  z i  is  even.  Since  B  is  odd,  y^  =  B  +  xi"zi 
is  odd  if  0  <  i  <  n^;  and  yn  is  even.  H^nce,  a±  *  1  for  0  <  i  < 
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and  a  *0. 
nl 


Q.E.  D. 


Theorem  3.2.2:  x^  >  y^  for  n  >  i  >  n^. 


k  mod  n2 

Proof:  Let  S'  »  S  +  2  .  Then  there  are  q  +  1  <  ng 

v. 

distinct  terms  in  S'.  It  is  easy  to  see  that  S'  <  2  -2,  and 

ni 

(2  -1)  S'  +  2r  <  B.  Thus,  by  (3.2.8) 


x±  -  (2  1-1)S*  +  2r 


-  yt  +  (2  -1)  *  2 


k  mod  n. 


(3.2.11) 


x^  >  y^  for  n  >  i  >  n^. 


Q.E.D. 


Property  5:  For  1  ^  0,  s  0  if  i  1e  a  multiple  of  n^  or  n£. 
Proof:  Case  1:  i  =  sn^  for  0  <  s  <  n2-l.  Since  k  *»  n-i  «  (^-s) 


n^,  r  »  0.  Equation  (3.2.8)  can  be  written  as: 


yt  -  (2  -1)[2 


(k-n^)  mod  02 


+...+20]  +  2° 


Therefore,  y^  is  even  and  »  0. 

Case  2:  i  ■  pn2  for  0  <  p  <  n^-1.  By  (3.2.9), 
n0  (k-n„)  mod  n.  n  n 

y±  *  (2  -1)[2  2  +. . .+2°]  +  2° 


Thus,  y^  is  even  and  ■  0. 


Q.E.D. 
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Property  6:  For  i  ^  0,  c^=l  and  a^*0  if  i  is  a  multiple  of  n2« 

Proof:  By  theorem  3.2.2,  we  have  x^-y^  >  0  if  i  is  a  multiple  of 
n2»  Thus  "  y^* 

By  equation  (3.2.4)  and  property  5,  z ^  is  odd  and  y.^  is  even. 
Therefore,  x^  is  odd  and  ■  1. 

Q.E .  D. 

Property  7:  If  i  is  not  a  multiple  of  n2»  then  a^  =  for  i  >  n^. 

Proof:  By  theorem  3.2.2,  we  have  <*  x^  -  >  0.  Since  z^  is 
even  if  i  is  not  a  multiple  of  x^  and  y^  must  be  either  both  even 
or  both  odd.  Thus  a^  ■  c^. 

Q.E.D. 

Fast  Algorithm  for  the  Binary  Form  of  A 

The  results  presented  above  enable  us  to  devise  a  simple 
algorithm  for  finding  the  binary  representation  form  of  A. 

9 

Let  us  denote  the  binary  form  of  A  by  the  n-tuple 


A  *  (Vl  Vz  • 


•  ai  V 


Then  a  ,  »  a  _ 
n-1  n-2 


an-n1-n2+l  "  °* 


n 


1  . 


Notice  that  the  binary  form  of  C  ■  2  *  A  is  just  the  binary  form  of 

A  shifted  nn  places  to  the  left,  i.e.  c..  *  a..  Let  us  denote  C  by 

1+n^  l 

C  "  (cn-l  cn-2  -  C1C0J 


By  properties  1  and  4, 


(a 


n. 


(0,  1,  1,  ...,  1,  1) 


Thus  (°2n  *  c2n  -i ■  •••*  ^  )  *  (0,  1,  1,  •••»  1»  1)  (3.2.12) 

Now  the  values  of  the  n. -tuple  (a.  ,  a„  -,...,  a  ,,)  can  be  determined 

1  r  2n^  2n^-l*  *  n^+l 

from  equation  (3.2.12)  and  properties  6,  7.  Next,  the  n-tuple 

(c3n  »  c3a  _lt  ...,  c2n  +1)  can  be  set  equal  to  (a2n  ,  ...,  an  +1> •  Again, 
111  11 
the  n^-tuple  (.a^  ,  •••,  a2n  can  determined  from  (c3n  ,  c3n  ^,...» 

c2u  +1^  properties  6  and  7.  This  process  can  be  repeated  until  all 

the  n-tuple  (ea_^,  an_2*  •••»  a^»  ag)  are  determined.  Therefore,  we  have 

the  following  algorithm: 

1.  Set  a,  *  1  for  0  <  i  <  n,-l,  a  “0  and  1*0. 

2.  If  (j+1)  n^  <  n-n^-n^,  set 

a.  c±+n^  »  &i  for  jnx+l  <.  i  <.  (j+1)  nx 

b.  a^  n  0  if  n^  divides  i 

a^  "  if  n2  does  not  divide  i 
for  (j+1)  nj+1  <  i  <  (j+2)  nx 

c.  j  *  j+1  and  repeat  step  2. 

3.  If  (j+1)  n^  i  n-n^nj,,  set  ai  «  0  for  n-n^-n2+l£i<n-l,  Stop 


The  algorithm  provides  a  fast  way  of  finding  the  binary  representation 
form  of  A,  even  for  very  large  integers.  For  example,  let  n^-4,  n2**9, 
A«(2^-l)/(2^-l)  (2^-1) .  Follow  the  above  algorithm,  the  binary  form  of 

A  is  readily  obtained  as; 

6  5  4  3  2  1 

c  oioooiooouoouoomomi 

A  000000Q009Qj^010p011001100nj£Ull 

7  6  5  4  3  2  1 


The  Complementary  Symmetrical  Property 


We  shall  show  that  the  binary  form  of  A  is  complementary  sym¬ 
metrical.  Specifically, 

V(nl+n2)-i+at  ‘  1  £or  • 

Let  us  consider  the  binary  forms  of  A  and  D.  Since, 
n,  (n.-l)n0  (n..-2)n„  n, 


1  v  1 
D  -  A(2  -1)*2 


2  v  1 
+2 


2  “2 
...  +2  +1 


D  **  1  000  ...  01  00  ...  010  ...  01  0  ...  0111  ...  1 

V _ /  s — v - - — /  ^  > 


n2  bits 


n2  bits 


n2-nl 


(3.2.13) 


where  x  denotes  -x.  Let  D*  be  the  integer  obtained  from  D  by  reversing  the 
order  of  the  sequences  in  (3.2.13),  i.e.. 


D*  »  1  1  ...  1  10... 01  00... 01  0...01  QI_^oi 


n.j-1  n2-ni 


n„ 


(3.2.14) 


Also  let  A*  be  the  i^reger  whose  binary  form  is  the  binary  form  of  A  in 

nl 

reverse  order.  Then  it  is  easy  to  check  that  D*  ■  A*  (1-2  ) .  Now,  let 

us  consider  the  integer  D-D*.  By  (3.2.13)  and  (3.2.14), 

D-D*  -  2  11. ..11  00... 01  ll. . .2  «  3  Q0...0ll  00... 0  00... 0 

n1  n1  n^  n2-ni  (ni"2)  n2 

-  3  00.... 03  00 . . . 9 

n^  (n1-l)n2-n1 

(n.— l)n«  (n. -1) n.-n.  (n-— l)n^— n-  n. 

Thus,  D-D*  »  3-2  -3-2  -  3-2  i(2  -1). 

nl  nl 

Since  D  ■  A(2  -1)  and  D*  ■  A*  (1-2  ) ,  we  have 
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A  +  A*  =  (D-D*)/  <2  X-l) 
n-n.-n 

=  3*2 

n-nr«2+i 


(3.2.14) 


Recall  that 


A  =  (1,  a  .  ••♦»  ai >  1) 

n“nl  2"1’’  1 


A  s  (if  fl.)  ••**)  8  ,  ,  1 ) 

v  ’  1*  *  n-ttj-n.-l  ' 


Starting  from  the  lower  order  bits,  by  (3.2.14)  we  have, 


or< in  general 


a.  4-  a  .4-1=0  mod  2 

1  n-n1-n2“l 


a.  4-  a  ,.4-leO  mod  2  for  0  <  i  <  n-n,-n„ 

i  n-n^-n^i  1  2 

Therefore,  a,  is  the  binary  complement  of  a  ,  ,  . 

i  r  n-n^-^-i,  and  we  have: 

Property  8:  The  binary  form  of  A  is  complementary  symmetrical,  i.e., 

a,  4-  a  ,  =  1  for  0  <  i  <  n-n.-nn  . 

i  n-n^n^i  1  2 

Property  8  enables  us  to  save  half  of  the  effort  finding  the  binary 

form  of  A  .  The  procedure  for  finding  the  binary  form  of  A  stops  at  the  place 

when  a^,  i  =  ^-(n-n^-n^  is  found.  The  high  order  bits  are  then  obtained  by 

reversing  the  order  of  the  complement  of  the  lower  order  bits.  The  comple- 

36 

mentary  symmetrical  property  can  be  seen  from  the  binary  form  of  A  =  (2  -1)/ 

4  9 

(2  -1)(2  -1)  in  the  previous  example. 


«wawBgi91^ 
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3.3  Construction  of  Multiple-Error-Correcting  AN-Codes 

>  •  ■ 

Let  A  be  the  generator  of  a  binary  cyclic  arithmetic  (AN)  code  of 
length  n.  Then  A  is  a  positive  integer  that  divides  2n-l.  Let  AB  =  2n-l,  for 
each  N,  0$  NS  B-l,  AN  is  a  code  word  of  the  code.  In  this  section,  we  shall 
consider  the  class  of  cyclic  AN-codes  whose  generator  A  is  of  the  same  form  as 
in  Section  3.2: 

A  *  ^n-l)/(2ni-l)(2n2-l)  (3.3.1) 

where  n^,  n^  are  relatively  prime  integers  and  n  =  nin2*  In  tlie  following 
discussion  we  shall  assume  that  >  n^ 

Any  code  word  AN  can  be  expressed  in  the  following  radix-2  or  binary 

form: 

AN  =  a  2n_1  +  an  2n"2  +  . . .+  a.  21  +  a. 
n-i  n-2  1  0 

0*0-1’  an_2>  •••  >  ai»  a0>  (3.3.2) 

where  a^  =  0  or  1  for  0<  i$  n-1.  From  the  definition  of  arithmetic  weight  of 
an  integer,  one  can  easily  derive  the  following  result. 

Lemma  3.3.1:  If  the  binary  sequence  of  AN  in  (3.3.2)  can  be  divided  into 
d  disjoint  subsequences  (each  subsequence  consists  of  a  certain  number  of  con¬ 
secutive  digits),  such  that  each  subsequence  contains  at  least  one  1  and  one  .  > 
then  the  arithmetic  weight  of  AN  is  at  least  equal  to  d. 

Given  A  and  B  such  that  AB  «  2n-l,  it  can  be  shown  [12]  that  the 
coefficients  a^s  in  (3.3.2)  are  given  by 


ai  5  (N*2  mod  B)  mod  2 


(3.3.3) 


Thus,  a  =  1  if  N2n  *  mod  B  is  odd  and  a.  =  0  if  N2n  *  mod  B  is  even. 

1  V  '  ”2 

For  the  code  considered  in  (3.3.1),  B  =  (2  -1)  (2  -1).  Thus-, 

n  n  nx 

2  (2  -1)  =2-1  mod  B 

in„  n.  n. 

2  (2  -1)  s  2  -1  mod  B,  for  any  positive  integer  i  (3.3.4) 

s  n2 

Lemma  3.3.2:  If  N  a  2  mod  (2  -1)  for  0  <  N  <  B-i  and  0  $  s  <  n2>  then 

the  arithmetic  weight  of  the  code  word  AN  generated  by  A  in  (3.3.1)  is  at  least 
equal  to  n^. 

Proof:  Let  a^'s  be  the  binary  coefficients  of  AN  in  (3.3.2).  By  (3.3.3) 

in  +n  -s 

a  ..  ,  .  5  (N2  4  mod  B)  mod  2 

n-(in2+n1-s) 

and 

in  -s 

a  /j  \  =  (N’2  mod  B)  mod  2 
n-(in2-s) 

where  0<  i<  n^-1,  and  the  subscripts  i  of  a^  are  taken  to  be  i  modulo  n. 

in„+n. -s  in  -s 

Let  x^  =  (N2  mod  B)  and  y  =  (N2  mod  B),  then  0  <  x^,' 

°1  n2  s  n2 

y±  <  B.  Since  B  =  (2  -1)  (2  -1)  and  N  s  2  mod  (2  -1), 

in_+n.-s  n_ 

xt  s  N2  mod  (2  -1) 


nl  n2 

5  2  1  mod  (2  -1) 


Similarly 
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in  -s  n 

yt  ■  N2  L  mod  (2  -1) 


s  1  mod  (2  -1) 


or 


Thus 


n„ 


yi  =  q2i  (2  1}  +  1 


n  n 

Xi  "  yi  “  (qli  “  q2i}  (2  _1)  +  (2  -1) 

n. 


“1  nl  no 

where  0  <  qu>  q u  <  2  -1,  for  0  <  x^  y±  <  B  =  (2  -1)(2  -1). 

hand, 


(a) 


On  the  other 


Thus 


or 


in  +n  -s  in  -s  in  -s  n 
N2  -  N2  L  =  N2  (2  -1) 

n 

=  N2~S(2  -1)  mod  B 

nl 

s  2  -1  mod  B 


Xi  “  yi  S  2  m0<*  ® 


n  n  n 

\  -  yt  =  q(2  -1) (2  Z-l)  +  (2  -1);  q  =  0,  or  -1 


by  (3.3.4) 


(b) 


nl  nl 

From  (a)  and  (b)  we  conclude  that  q(2  -1)  =  qljL  -  q2±.  But  jq^  -  q9J  <  2  -1, 


*li  M2i 1 


q  must  be  0.  Thus  x^  -  y^  =  2  -1.  This  implies  either  x^  =  odd  and  y^  =  even 
or  x^  =  even  and  y^  =  odd.  In  othei  words,  for  each  i,  one  and  only  one  of  the 
pair  o£  binary  digits  an-(ln2-s)  is  e’ual  t0  1  'or  °>- 

Now,  the  sequence  of  n  digits  of  AN  in  (3.3.2)  can  be  divided  into  n^ 
subsequences,  each  subsequence  contains  consecutive  digits.  Therefore;  the 
digits  an„^n  +n  _g)  an^  an-(in  -s)  are  *n  same  subsaquence  for  a  fixed  i. 
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Each  subsequence,  then,  contains  at  least  one  1  and  one  0.  By  Lemma  1  the 

arithmetic  weight  of  AN  is  at  least  n. . 

1  Q.E.D. 

n9 

Lemma  3.3.3:  If  N  h  -2  mod  (2  -1)  f or  0  <  N  ^  B-l  and  0  <  s  <  n2>  the 
arithmetic  weight  of  the  code  word  AN  generated  by  A  in  (3.3.1)  is  at  least 
equal  to  n^. 

Proof:  This  lemma  can  be  proved  in  a  similar  way  as  in  Lemma  3.3.2. 

Q.E.D. 

n0  s  n2 

Lemma  3.3.4:  If  N  #  0  mod  (2  -1)  and  N  4  ±  2  mod  (2  -1)  for  0  <  N  <  B-l 
the  arithmetic  weight  of  AN  generated  by  A  in  (3.3.1)  is  at  least  equal  to  n^. 

Proof:  Let  V  =  AN,  C  =  V(2  -1)  and  W(x)  denote  the  arithmetic  weight  of  x 

Since  the  arithmetic  weight  of  integers  satisfies  triangular  inequality, 

n.  n 

W(C)  =  W[V(2  -1)]  2s  W(V2  l)  +  W(V)  =  2W(V) 

W(C)  <  2W  (V) 

n„  (n. -l)n„  (n.-2)n2 

Let  Nj  =  N  mod  (2  -1),  s  C  mod  (2-1).  Then,  Cj  =  Nj2  +^2 

n9 

+  ...  +  N^2  +  N^.  Thus  the  binary  form  of  consists  of  n^  replicas  of  the 

binary  form  of  N^.  Since  is  nonzero  and  N^  f  +  2s,  the  binary  form  of  N^ 

contains  at  least  two  l's  and  two  0’s.  Therefore,  the  modular  weight  [36]  of 

C. ,  W  (C.)  2;  2n, .  Hence, 
lml  1  •  ’ 

2W(V)  *  W(C)  2:  W  (C.)  a  2n. 

mi  I 

W(V)  nx 


Q.E.D. 


Thaorem  3.3.1:  The  minimum  arithmetic  distance  of  the  code  generated  by  A 
in  (2.3.1)  is  equal  to  n^  for  n^  <  n^. 

Proof:  Let  us  assume  first  that  N  s  0  mod  tf”2-!),  0  <:  N  £  B-l.  Then 
N  =  q<2  2-l),  0<q<:  2\l,  and 


nln-j 
2  1  Z-1 


2  -1 


(n?-l)n  (n„-2)n. 

»q2  2  1  +  q  2  2  2+. 


•  +  q’2  +  q 


By  lemma  3.3.1,  W(AN)  ^  «2  >  nL  . 

For  other  nonzero  values  of  N,  lemmas  3.3.2,  3.3.3,  and.  3.3.4  assure 
that  the  code  word  AN  he.  en  arithmetic  weight  of  ^  or  greeter.  Therefore,  th 
arithmetic  distance  of  th.  code  i»  at  leaat  equal  to  v  However,  the  weight  o 

the  code  word  A (2  2-l)  is  exactly  equal  to  ly  Thu.  the  arithmetic  distance  o 
the  code  is  equal  to  n^. 

Q.E.D. 

The  code  we  constructed  above  is  a  multiple-error-correcting  cyclic 
AN-code.  The  synthesis  of  this  class  of  codes  is  straightforward,  we  merely 
uhcose  two  proper  relatively  prime  integers  ^  and  n^  the  integer  tf"1"2-!)/ 

(2  -1)(2  -1)  generates  a  code  with  minimum  distance 

3.4  Root-Distance  Relation  of  p1?  -Codes 

If  the  block  length  n  is  a  product  of  two  distinct  primes  ?1  and  p  , 

the  code  is  called  a  p^-code.  \  p^-code  is  an  ideal  in  the  ring  of 
integers  modulo  2  1  2-l 
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pLpo 

Since  the  polynomial  x  -1  can  be  factored  into  cyclotomic  polynomials 
over  the  rational  field  as: 

pipo  1  Pi  *o  pipo 

x  1  -1  =  Q  (x)Q  L(x)Q  Z(x)Q  Z(x)  (3.4.1) 

We  can  substitute  2  for  x  in  (3.4.1)  and  get 

pTpo  i  pi  po  pipo 

2  1  -1  =  Q  (2)Q  i(2)Q  Z(2)Q  1 (2)  (3.4.2) 

where  Q^(2)'s  are  called  the  cyclotomic  factors  of  2n-l. 

We  define  a  "cyclotomic  generator"  of  arithmetic  codes  as: 


Definition  3.4.1:  A  cyclotomic  generator  of  an  arithmetic  aode  is  a 
product  of  some  cyclotomic  factors  of  2n-l. 

By  this  definition,  the  generator  of  the  codes  considered  in  Section 
3.3  is  a  cyclotomic  generator. 


The  only  possible  cyclotomic  generators  of  a  p^p^-code  are: 


1  Pi  ^2  ^1  ^2  pi  pip2 

Aj_  =  Q  (2),  A2  -  Q  L(2),  A3  =  Q  t(2),A4  =  Q  X(2)Q  Z(2),A5  -  Q  L(2)Q  Z(2), 

P2  P1P2  pip2  pip2 

Afi  -  Q  A(2)Q  i2) ,  A.7  =  Q  1  ^(2)  and  Ag  =  2  -1.  Both  and  Ag  are  trivial 

generators. 

P1  P1  P1 

Since  Q  (2)  =  (2  -l)/(2-l)  =2  -1,  the  minimum  distance  of  the 

P1 

code  generated  by  A^  *  Q  (2)  is  2.  Similarly,  the  minimum  distance  of  the  code 

generated  by  A.  is  also  2.  These,  are  the  single  error-detecting  codes.  For 
P1  p2 

A  =  A^  =  (2  -1)(2  -1),  the  minimum  distance  of  the  code  is  given  by  the 


following  theorem: 
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Theorem  3.4.1:  The  minimum  distance  of  p^-code  generated  by  A4  is 

(a)  3  if  pxp2  =  6; 

(b)  4  if  Plp2  £  6. 


Proof;  With  no  loss  of  generality  we  may  assume  P2  >  Pr  hence  p2  i  3. 
Furthermore,  PjP2  3  6  if  and  only  if  =  2  and  ?2  =  3.  It  is  obvious  that  A 4 

cannot  divide  21  hence  there  .re  no  codevord.  of  weight  1.  Codeword,  of  weight 
2  will  be  of  the  fora,  2  V  ±  i) ,  where  m  <  p^.  For  ^  ,  2k(2m-l)  h.ve 

(2  l-l)|(2”-D,  hence  pjm.  Similarly,  p2|„  hence  p^m.  But  p^  i.  gre.ter 
th.n  .  contradiction.  If  y  -  2k(2\l>,  then  we  may  write  a,  =  ^  +  ^ 

«  -l)|cz“+l)  then  implies  2  ’-llz^+l,  ,n  impossibility.  Hence  .11  nonzero 

codewords  of  the  code  generated  by  A  =  o^2  o1  ,P1. , 

B  ueu  uy  »  -2  -2  +1  are  at  least  of  weight  3. 


,4. -2 


For  the  case  Pjp2  *  6,  p1  *  2  and  p2  »  3.  Hence  A^  ■  25-23-22  +  1 


£  * 

-2+2  +1.  Consequently  A4  generates  a  code  of  minimum  distance  3 


p  F°p  the  Ca8e  P1P2  *  6>  iC  that  p2  *  Pl  +  2,  hence  A4  « 

2  -2  +1  is  of  weight  4.  The  codewords  of  weight  3  will  take  the  form 


2  (2  +2  +1) ,  where  m2  £  m^-2 . 

Writing: 

“i  *  hi^n 

where  0$  r^  <  Pl 

"l  *  <*12p2+r12 

where  0  <  r12  <  p2 

"2  *  ’21pl+r21 

where  Os  r21  <  p 

“2  *  ’22|,2+r22 

where  0  <  r22  <  p2 

we  may  deduce  that 

pi  ,  ml 

2  l-l  | 2  +  2  2  + 

„  pl  rii 

1  implies  2  -1)2  LL+  2 
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and 


Pn  ®i  Po  ,  rio  r22 

2  -lj2  L+  2  Z  +  1  implies  2  -l|2  +2+1 


k  ml  m2 

*•  A  t 


The  only  possibilities  are:  A^N  =  2^(2  *-2  ‘-l)  where  r^  =  ^2  = 
•12  =  rll 


k  ml  m9 

r.,  -  ru  =  1  and  A.N'  =  2  (2  +2  -1}  where  r21  =  rn  =  pj-1;  r22  -  r12  =  p2~l. 


m  ® 

When  A4N  -  2k(2  1-2  2-l),  we  have: 

k.t,p!p2+l  _28'P1P2  ml) 


A4N  -  2  (2 


This  is  clearly  impossible  as  s  Ss  1  and  m2  <  m^-2  would  imply  t  i  2.  When 


k  ml  rao 

A4N'  -  2  (2  +2  -1),  we  have: 


A4N*  =  2  (2 


k__t.PlP2-l  8-P,P,-l 


+  2  «  -1) 


Again  si  1  and  m2  £  n^-2  would  imply  t  i  2,  a  contradiction.  This  completes 
the  proof  of  theorem  3.4.1. 


Q.E.D. 


P1  P1P2 

For  A  -  A5  »  Q  (2)Q  *(2),  we  have: 

P1P2  ,_P2. 


PiP 


1**2 


A5  =  (2  -1)/Q  (2)  =  (2  -1)/ (2  -1) 


<P1-!)p2  (P1-2)P2 


+2 


+  . .  ♦  +  2  +1 


As  indicated  by  Erosh  [23],  the  minimum  distance  of  a  code  generated  by  an 

klk2  kl 

integer  of  the  form  (2  -l)/(2  -1)  is  k,,.  Therefore,  the  p^Pj-code 

generated  by  A,.  is  of  distance  p^.  Similarly,  the  code  generated  by  A^  is  of 


distance  p,,. 
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PiP,  P1P0  Pi  P? 

For  A  »  A7  =  Q  1  Z(2)  =  (2  -l)/(2  -1) (2  -1),  the  code  is  a 

/ 

special  case  of  the  code  mentioned  in  Section  3.3.  Since  (p^^)  =  ^  an<* 

P1  <  P2’  by  theorem  3.3.1,  the  minimum  distance  is  p^. 

From  the  above  results,  some  interesting  properties  about  the  root- 

distance  relation  of  arithmetic  codes  can  be  brought  to  light. 

Let  us  consider  the  polynomial  xn-l  over  the  rational  field.  If  u 

is  a  primitive  n-th  root  of  unity  [48],  then  all  roots  of  x  -1  are  powers  df  u 

n  . 

and  we  have  xn-l  =  Also,  x  -1  can  be  factorized  into  cyclotomic 

polynominals  over  the  rational  field  as: 

xn-l  =  tt  Q^(x)  (3.4.3) 

d/n 

Hence,  the  roots  contained  in  each  Qd(x)  are  powers  of  u.  Due  to  the  property 
of  cyclotomic  polynomials  [6],  the  roots  in  Qd(x)  are  the  u^*s  where  i's  satisfy: 
(n,i)  =  n/d. 

Substituting  2  for  x  in  eq.  (3.4.3),  we  get  2n-l  =  n  Q  (2)«  Let  A 

d/n 

be  a  cyclotomic  generator,  then  A  equals  a  product  of  some  cyclotomic  factors, 
say: 

dl  d2  d3  ds 

A  =  Q  \2)Q  Z(2)  ■i(2)....Q  S(2) 

We  define  the  corresponding  polynomial  of  the  generator  A  as: 

di  d9  d 

A(x)  =  Q  (x)Q  (x) - Q  S(x) 

Then  the  roots  of  the  integral  generator  A  can  be  defined  as  follows: 
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Definition  3.4.2:  The  roots  of  a  cyclotomic  generator  A  are  defined  as  the 
roots  of  its  corresponding  polynomial  A(x) . 

This  definition  of  the  roots  of  A  was  first  proposed  by  Chien  and 
Hong  [15].  Since  all  the  roots  of  A  are  powers  of  u,  the  number  of  consecutive 
roots  in  A  is  defined  as  the  largest  number  of  consecutive  powers  of  u  contained 
in  A(x) . 

For  the  p^p^-code,  the  root-distance  relation  can  be  easily  obtained. 

PiPo  i  Pi  Po  PiP? 

The  cyclotomic  factors  of  2  -1  are  Q  (2),  Q  (2),  Q  L (2)  and  Q  (2).  If  u 

is  a  pjp^-th  root  of  unity,  then  with  a  little  calculation,  we  can  summarize  the 

root-distance  relation  as  follows: 


Generator 

No.  of  consecutive  roots 

Minimum  distance 

P1 

A2  =  Q  L(2) 

1 

2 

p2 

A3  =  Q  l2) 

1 

2 

P1  p2 

A4  -  Q  i(2)Q  Z<2) 

*if  p^  **  2,  p2  f  3 

3 

4 

*if  P^»P2  are  °di 

2 

4 

*if  PXP2  =  6 

3 

3 

P1  pip2 

A5  -  Q  i(2)Q  L(X) 

pr1 

P1 

P2  P1P2 

A6  -  Q  Z(2)Q  Z(2) 

V1 

P2 

P1P2 

A?  -  Q  Z(2) 

pri 

P1 

From  the  above 

tabulated  results,  we  see  that 

for  all  possible 

cyclotomic  generators  of 

a  pjf>2-code  (except  the  case  p 

^P2  =  6),  the  minimum 

distance  of  the  code  is  at  least  one  greater  than  the  number  of  consecutive  roots 
contained  in  the  generator — the  same  root-distance  relations  as  the  BCH  theorem 
for  polynomial  codes.  As  for  the  case  of  P^2  =  6,  we  have  the  following  argument. 

Since  26-l  =  Q1(2)Q2(2)Q3(2)Q6(2) ,  where  Q1(2)  =  1,  Q2(2)  =  3,  Q3(2)  =  7 

g 

and  Q  (2)  ■  3,  the  generator  is  21.  The  code  generated  by  A^  is  of  distance 

2  3 

3,  which  equals  the  number  of  consecutive  roots  contained  in  Q  (x)Q  (x) .  How- 

3 

ever,  the  generator  A^  =  3.7  can  also  be  considered  as  the  product  of  Q  (2)  and 
Q  (2)  since  Q°(2)  =  Q  (2)  =  3.  The  number  of  consecutive  roots  in  QJ(x)Q  (x)  is 
only  2,  which  is  indeed  one  less  than  the  minimum  distance  of  the  code.  The  situa¬ 
tion  that  a  cyclotomic  generator  can  be  expressed  as  different  products  of  cyclo- 
tomic  factors  can  only  happen  when  p.^  =  since  Dickson  [21]  has  shown  that 
Qi(2)  f  Q^(2)  for  all  i  f  j,  except  i  =  2  and  j  =  6. 

The  analogies  between  cyclic  parity-check  codes  and  cyclic  AN-codes 
have  been  investigated  by  many  authors.  The  root-distance  relationship  discussed 
above  indicates  another  analogy.  It  is  our  hope  that  the  closeness  of  all  these 
analogies  would  lead  us  to  the  discovery  of  a  class  of  cyclic  AN-codes  which  is 
analogous  to  the  large  and  powerful  Bose-Chaudhuri-Hocquenghem  codes.  Then,  for 
all  practical  purposes,  the  synthesis  of  multiple-error-correcting  code  can  be 
achieved  through  the  root  conditions  of  the  code  generator. 


3.5  Bounds  on  Error  Control  Capability  of  Cyclic  AN-Codes 


Since  the  generator  A  of  cyclic  AN-code  of  length  n  is  a  proper  divisor 
of  2n-l,  we  can  write  2n-l  =  AB,  with  B  >  1.  As  briefly  mentioned  in  the  intro¬ 


ductory  section  of  this  chapter,  the  well-known  "middle-third"  region  of  B  is 
defined  as  the  set  of  consecutive  integers  x,  with  B^  3  x<  2  B.  If  we  denote 


v»»?*73V'  w-'W?*1 
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the  "middle-third"  region  of  B  by  Mg,  then  there  are  !~(B+1)/31  integers  in  Mg, 
where  fx]  denotes  the  largest  integral  value  of  x. 

It  has  been  shown  [12,17]  that  the  arithmetic  weight  of  any  codeword  is 
related  to  the  integers  in  the  region  Mg.  This  result  can  be  briefly  summarized 
in  the  following: 

Theorem:  The  arithmetic  weight  of  a  codeword  AN  is  given  by  the  number 
of  residues  c  Mg,  where,  Xt  *  (N‘21)  mod  B,  for  i  =  0,1,2,...,  n-1. 

From  the  above  theorem,  it  is  seen  that  the  minimum  distance  of  a  code 
can  be  evaluated  by  counting  the  number  of  residues  in  the  region  Mg. 

To  calculate  the  distance,  we  first  divide  the  codewords  into  several 
disjoint  orbits  [17],  for  each  orbit  we  count  the  number  of  residues  in  Mg,  and 
get  the  arithmetic  weight  of  each  orbit.  Then  the  distance  equals  the  minimum 
weight.  Clearly,  this  evaluation  will  be  quite  complex  for  large,  and  composite 
B's.  Also,  since  we  do  not  have  a  closed-form  formula  for  the  distance  of  a  code, 
the  above  result  does  not  suggest  a  synthesis  procedure  for  finding  any  multiple- 
error  correcting  arithmetic  codes,  we  can  only  calculate,  through  a  computer- 
programmed  procedure,  the  error-correcting  capability  of  a  given  code.  There¬ 
fore,  the  result  is  a  theoretical  analysis  of  the  structure  of  the  arithmetic 
codes  rather  than  a  practical  formulation. 

In  this  section,  we  shall  apply  some  number  theoretic  concepts  to  ob¬ 
tain  the  upper  bound  and  the  lower  bound  on  the  error-control  capability  of  any 
cyclic  AN-codes.  These  bounds  suggest  the  criteria  for  us  to  choose  proper 
parameters  in  constructing  an  AN-code. 
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3.5.1  Lower  Bound  on  the  Distance 

Let  the  B  integers  { 0 ,  1,  2,  ...,  B-l}  be  divided  into  three  regions 
Lg>  Mg,  and  Ug,  where: 

Lg  =  "lower  third"  region  =  fx|0  <  x  <  f^l  } 

Mg  =  "middle  third"  region  =  {x  |  <  x  <  F—l  } 

o 

Ug  =  "upper  third"  region  =  [x|  IjBl  <  x  <  B} 

Let  us  define  the  following: 

Definition  3.5.1:  For  any  integer  N,  0  <  N  <  B,  the  smallest,  non-negative 
integer  i  such  that  (2^-N)  mod  B  e  M„  is  defined  as  the  index  of  N»  and  the  inte- 

D  — — 

ger  N  is  said  to  have  index  i. 

Suppose  that  AN  is  a  codeword,  then  0<  N  <  B.  it  ft-jM  ,  then  the 

smallest  power  of  2  that  yields  2^N  mod  BeMg  is  clearly  zero,  hence  N  has  index 

0.  For  NeLg  or  NeUg,  the  index  is  certainly  not  zero.  Let  us  firsc  consider  the 

integers  in  Lg.  Except  the  integer  0,  such  integer  in  Lg  has  a  nonzero  index. 

Denoting  the  index  of  xeLg  by  i^,  we  have  the  following: 

Lemma  3.5.1:  The  index  i,  =  ["log-  . 

1  2  3 

*1  *1~^ 

Proof:  Since  2  *1  mod  BeMg  and  2  ’1  mod  BgLg,  we  have: 

V1  B  il 

2  1  <  r|i  <  2  1 

or: 

i^l  <  log2  r|l  <  ^ 

thus: 


il  ri°g2  3  ^  ’ 


Q .E.D. 
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Lama  3.5.2:  For  xiyeLB,  ^  iy  if  !  ,  >  j  >  „.  Purthcraore,  if,:!,, 

then  i  >  i  . 
x  y 

i  Pr00f:  Since  Xiy«v  by  the  definition  of  index,  the  integers  2lx*x  and 
2  y  are  less  than  B,  thus  we  have; 


2  x  -xeMB  and  2  y  -yeMg  . 

Also,  by  the  definitions  of  LQ,  ^  and  Ug,  we  have: 

r-B  1  ^  2  x  x  =  2  x  •  (2x)  <^4  1 


(3.5.1) 


r^Bl  <  2  y  -y  r|Bl 
Hence,  if  2  x  >  y  >  x,  we  have: 


(3.5.2) 


A  i-w  i  -1  i  -1 

2  *y  >  2  “'•x  =  2  X  (2x)>2X  -y 


(3.5.3) 


From  (3.5.1),  (3.5.2)  and  (3.5.3),  it  is  seen  that  the  integer  2  x  'y  is  either 


in  Lb  or  in  Mg. 


*y  eV  then  \r  "  Vi?  ^ile  if  2  x  -ycL  ,  then 

J  D 


2  ’yeMg  and  i  =  i  ,  hence  i  *  i  . 

x  y  x  y 

If  y  i  2x,  then  we  have: 


i  i  +1 

i  y*y  2  2  y  *x 

.  i  +1  i  +1 

hence  either  2  ^xeMg  or  2  y  -xeL  . 

i  +1  ±  +1 

If  2  y  -xeM .  then  i  +1  =  i  ,  while  if  2  y  • 


hence  i  >  i  , 
x  y 


y  xs 


■xgL  then  i  +2  <  i  ; 

n  w  v* 


y  x’ 


Q.E.D. 
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n  - 1 

Lemma  3.5.3:  The  residue  2  mod  B  is  in  if  B  f  1. 
Proof:  Let  2n  *  a  k  (mod  B),  then 


Since 


2n_1  =  m  •  B  +  k 


for  some  m. 


2n-l  =  AB,  both  A  and  B  are  odd,  and 

2n  =  AB  +  1 


hence, 

_n-l  AB+1 
2 


»H  .  (A-l)B  _  1  _ 

=  -j-  +  — 2*~  ~  j  (B+l)  mod  B 

But 


-  |  =  ^  [3B+3-2B]  =|[B+3], 

and 

'f  -  -X  -  |  [4B-3B-3]  =  \  [B-3J 
Hence,  for  Bi  3,  2  mod  B  =  — (B+l)  cM  . 

Z  B 


(3.5.4) 

i 

(3.5.5) 


{ 


Q.E.D. 


Since  the  integers  in  U„  are  congruent  to  the  negatives  of  the  integers 

B 

in.  L  ,  we  can  use  the  same,  argument  as  in  lemmas  3.5.1  and  3.5.2  to  obtain: 

B 


Lemma  3.5.4:  For  xgU-,  the  index  i  equals  to  index  i  ,  where  B-x$L  . 

B  X  B-X  B 

Therefore,  in  the  range  ft),  B-ll  ,  the  integers  1  and  B-l  have  the 
largest  index  Tlog^  —1  .  With  the  above  lemmas,  we  now  derive  a  lower  bound  on 
the  minimum  distance  of  the  cyclic  AN-codes. 


rr^  7K^r'Ti«TVT^"--’'  'r'T  " 


-'figm I '•',*-’TKT7VZ~!  ^*TT*C<**' \*IX 
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Let  us  consider  an  orbit  of  n  residues: 

(N  mod  B),(2N  mod  B) ,  (22N  mod  B) , .  (2n_1N  mod  B)  (3.5. b) 

Suppose  that  the  residue  (2^N  mod  B)eMg.  By  lemma  3.5.1,  there  exists 
another  residue  (2k+mN  mod  B)eMg,  with  l<m<i^+l,  where  i^=*  og . 

The  rest  of  the  residues  are: 

(2k+ntflN  mod  B),(2k+nrf'2N  mod  B)  , . . . ,  (2n-1N  mod  B),(N  mod  B), 

(2k_;lN  mod  B)  (3.5.7) 


By  lemmas  3.5.2  and  3.5.3,  we  see  that  there  exists  at  least  one  residue 
for  every  i^+1  consecutive  residues  in  the  set  (3.5.7). 
Since  the  arithmetic  weight  of  the  code  word  AN  equals  the  number  of  resi¬ 
dues  of  the  orbit  (3.5.6)  which  are  in  Mg,  we  have  proved  the  following 
theorem: 

Theorem  3.5.1:  The  minimum  distance,  D^,  of  a  cyclic  AN-code  of 
length  n  is  at  least  as  large  as 


which  is  in 


Dm-  r^r^r^-l  +  2  (3-5*8> 

The  tightness  of  the  bound  in  theorem  3.5.1  can  be  improved  if 
we  consider  the  following: 

If  B*3,  then  i^=0.  In  this  case  there  are  only  two  codewords  in 
the  code,  the  minimum  distance  of  the  code  is  trivially  seen  to  be  n/2.  In 
the  following,  we  consider  the  case  that  B?*3: 

Lemma  3.5.4:  Let  the  exponent  of  2  modulo  B  be  denoted  by  e(B).  If 
e(B)“n,  and  if  1  and  B-l  are  not  contained  in  the  same  residue  set 
(2^  mod  B;  0<i<n},  then 


D 

m 


■2 
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Proof:  Since  1  and  B-l  are  in  distinct  residue'  sets,  for  every  i^ 
consecutive  residues  in  the  residue  set  (2)  there  must  exist  at  least  one 
residue  in  Mg. 

Lemma  3.5.5:  If  e^J^n  and  2  =-l  mod  B  for  some  0<k<n,  then  n^k. 
Proof :  We  have 


e(B)  =  1 

(mod  B) 

k  5  -1 

(mod  B) 

2k  =  1 

(mod  B) 

By  definition  of  e(B)  ,  n  is  a  divisor  of  2k,  If  n  is  odd,  then  n|k;  this 
is  impossible  since  k<n.  If  n  is  even,  then  (n/2)|k  or  k=t*(n/2)  for 
some  integer  t.  Since  0<k<n,  t=l.  Therefore,  n=2k. 

Lemma  3.5.5  states  that  if  1  and  E-l  are  contained  in  the  same 
residue  set  {2  mod  B;  CKL<n},  the  length  of  the  code  must  be  even.  Any 
residue  set  {N'2'1'  mod  B;  0<i<n}  can  then  be  divided  into  two  subsets, 

SL  =  IN  mod  B,  N*2  mod  B,  . . . .  N.2(n/2)_1  mod  B) 


S2  =  {N,2^n/2'  mod  B 


,  N-211"1  mod  B} 


The  elements  in  d™  tne  additive  inverses  in  the  ring  of  integers  modulo 
B  of  the  elements  in  S^.  These  two  sets  have  the  same  number  of  residues 
that  are  in  the  region  Mg.  Therefore,  we  have 

Lemma  3.5.6:  If  e(B)=n=even,  and  2^n^=-l  mod  B,  then 
<n/2)-i  -2 

D  >  - — — •  +2  x2 

m  —  I  i^  * 

it  should  be  mentioned  that  (n/2)-i^-2  is  always  non-negative  as  long  as 
the  codes  have  error  correcting  capability.  This  is  true  because  that 
(n/2)-i^\L  ana  the  equality  hold.i  only  when  the  code  is  of  distance  2. 


If  the  exponent  of  2  modulo  B  is  not  equal  to  n,  then.  e(B)  is 
a  factor  of  n,i,e.,  n=>e(B)*q.  In  this  case  any  residue  set  {N*2^  mod  B, 

0^i<n}  can  be  divided  into  q  identical  subsets. 

Lemma  3.5.7:  If  n=e(B)-q,  1  and  B-l  are  not  contained  in  the  same 
residue  set  {21  mod  B,  CKL<e(B)},  then 

e(B)-i  -2 

(a)  | - ^ —  +2  j  xq  for  eCBj-i^l 

(b)  =  q  for  e(B)-i^=l 

Proof: 

(a)  Since  e(B)-i.jf\L,  we  have  e(B)-i^-2>^-l.  If  e(B)-l^>l,  then  e(B)-i^-2^0. 
From  lemma  3.5.4  and  the  fact  that  any  residue  set  contains  q  identical 
subsets,  part  (a)  of  the  lemma  follows  clearly. 

(b)  If  e(B)-i^=l,  then  e(B)-i^-2=-l.  This  implies  that  among  the  e(B) 
residues  2N  mod  B,  2^N  mod  B,...,  2e ^  mod  B,  there  is  only  one  residue, 
namely,  2e^  mod  B,  which  is  in  Mg.  The  NAF  of  the  generator  A  is  seen 
to  he 

e(B)  bits  e(B)  bits  e(B)  bits 

00. . .  .01,00.  t  .01, . 00..  .01 

and  the  minimum  distance  of  the  coda  is  q.  Q.E.D. 

If  n=e(B)*q,  and  if  1,  B-l  are  Doth  contained  in  the  same  resi¬ 
due  set  {2*  mod  B,  0^i<e(B)},  t'ren  by  a  similar  argument  as  in  lemmas  3.5.5 
and  3.5.6  we  have: 

Lemma  3.5.8:  If  n=e(B)’q  and  2e^^=-l  mod  B,  then 
e(B)/2-i,-2 

D»i  [ - r^-+2l-2., 

Again,  by  using  che  same  argument  as  in  lemma  3.5.7,  we  have 


Lemma  3.5.9:  If  n*e(B).q,  2e^^=-l  mod  B  and  e(B)/2-ij“l,  then 


D  -  2q 
m  H 


In  light  of  the  above  lemmas,  the  determination  of  the  lower 
bound  on  the  minimum  distance  of  any  cyclic  AN-code  can  be  summarized  as 


follows : 


Theorem  3.5.2:  For  a  cyclic  AN-code  of  length  n. 


(a)  n=e(B) 


Dm  1  {r«n/2)-i1-2)/i1l+2}.2 

if  n«even  and  2^n  2^s~l  mod  B 

\  1  f(n-ir2)/iil  +2 


otherwise 


(b)  B-3 


D  -  n/2 
m 


(c)  n«e(B)*q 


\  >  {  f(e (B) / 2-i.j-2) /i^*|  +2}-2q 

;if  e(B'«even  and  2e^B^2=-l  mod  B 
**  ^m  °  ^q  ;if  e(B)*even,  2e^B^^=-l  mod  B, 


and  e(B)/2-lj»l 


Dffl  >.  {  f(e(B)-i1-2)/i1'J  +2}.q 

;if  2e(B)/2*-l  mod  B 


**  D  *  q 
m  n 


;if  2e(B)/2^-l  mod  B  and  eCBj-ij-l 


For  any  given  cyclic  AN-code,  the  above  theorem  can  be  used  to 
determine  the  lower  bound  on  the  error  correcting  capability  of  the  code. 

The  calculation  is  very  simple.  We  present  in  the  following  the  lower 
bounds  (calculated  by  the  CDC  1604  computer)  of  all  possible  cyclic  AN-codes 
of  length  up  to  36: 
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I 

I 

I 

I 

I 

1 

I 

I 


n(code  length) 

B 

e(B) 

8 

3 

2 

8 

5 

4 

9 

7 

3 

10 

11 

10 

10 

2 

12 

7-13 

12 

12 

3.3*13 

12 

12 

3*3*5 

12 

12 

5*7 

12 

U 

13 

12 

12 

3*3 

6 

14 

43 

14 

14 

3 

2 

15 

7*31 

15 

15 

151 

15 

15 

7 

3 

16 

5*17 

8 

16 

51 

8 

16 

3*5 

4 

16 

5 

4 

16 

3 

2 

18 

3*3*7 

6 

18 

3*3*73 

18 

18 

7 

3 

18 

3*3 

6 

18 

7*3*19 

18 

18 

73 

9 

20 

3*11*5*41 

20 

20 

3*31*5*41 

20 

20 

5*3*31 

20 

20 

5*3*11 

20 

20 

5*11*31 

20 

20 

11 

10 

20 

3*5 

4 

20 

5*11 

20 

20 

5 

4 

20 

3 

2 

21 

7*127 

21 

21 

7*337 

21 

21 

7 

3 

22 

683 

22 

22 

3 

2 

24 

3*3*5*7*17 

24 

24 

3* 3* 5* 7* -41 

24 

24 

5*7*31*241 

24 

24 

3*  3*5*13*241 

24 

24 

3*3*13*7*17 

24 

i,  D  (lower  bound)  D  (actual) 
i  m  m 


0  4  4 

14  4 

2  3  3 

2  4  4 

0  5  5 

5  3  3 

6  3  3 

4  3  3 

4  3  4 

3  4  4 

2  4  4 

4  4  4 

0  7  7 

7  3  3 

6  3  3 

2  5  5 

5  4  4 

5  4  4 

3  4  4 

18  8 

0  8  8 

5  3  3 

8  3  4 

2  6  6 

2  6  6 

8  3  4 

5  4  4 

12  3  3 

14  3  3 

8  3  4 

6  4  4 

10  3  4 

2  8  8 

3  5  5 

5  4  4 

1  10  10 

f  10  10 

9  3  3 

10  3  4 

2  7  7 

8  4  4 

0  11  11 

11  3  3 

15  3  3 

16  3  3 

12  3  3 

12  3  4 


n(code  length) 

B 

e(B> 

24 

3*3*5*  7 

12 

24 

3*3*5*13 

12 

24 

3*3*13*17 

24 

24 

7*13*17 

24 

24 

5*7*17 

24 

24 

3*3*13 

12 

24 

3*3*5 

12 

24 

5*7 

12 

24 

13 

12 

24 

9 

6 

24 

5 

4 

24 

3 

2 

25 

31 

5 

26 

2731 

26 

27 

73 

9 

27 

7 

3 

28 

3*5*43 

28 

28 

3*5*127 

28 

28 

3*29*113 

28 

28 

3*127.29*113 

28 

28 

3*43*29*113 

28 

28 

3*29*113 

28 

28 

3*5 

4 

28 

5 

4 

28 

3 

2 

30 

3*3*11*7*31 

30 

30 

7*151*11*331 

30 

30 

9*31*151*331 

30 

30 

9*11*31*151 

30 

30 

7*11*31*331 

30 

30 

7*9*31*331 

30 

30 

7*9*11*151 

30 

30 

9*31*33.1 

30 

30 

9*31*151 

30 

30 

9*11 

30 

30 

9*7*11 

30 

30 

7*9*31 

30 

30 

9*31 

30 

30 

9*11 

30 

30 

7*11*31 

30 

30 

3*3 

6 

30 

11 

10 

30 

7 

3 

30 

3 

2 

32 

3*5*17 

8 

32 

3*5*257 

16 

32 

3*17*257 

16 

32 

5*17*257 

16 

Dm( lower  bound)  Dffl(actual) 


7  4  4 

8  4  4 

10  3  4 

9  3  5 

8  3  4 

6  4  4 

4  6  6 

4  6  8 

3  8  8 

2  8  8 

1  12  12 

0  12  12 

4  5  5 

10  4  4 

5  6  6 

2  9  9 

8  4  4 

10  4  4 

12  3  4 

18  3  3 

17  3  3 

12  3  4 

3  7  7 

1  14  14 

0  14  14 

13  3  3 

21  3  3 

22  3  3 

18  3  3 

19  3  4 

18  3  4 

14  4  4 

12  3  4 

14  3  4 

6  6  6 

8  4  5 

10  4  5 

7  5  6 

6  6  6 

10  4  5 

2  10  10 

2  12  12 

2  10  10 

0  15  15 

7  4  4 

11  4  4 

13  4  4 

13  4  4 
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n(code  length) 


32 

32 

32 

32 

32 

33 
33 

33 

34 

34 

35 
35 
35 

35 

36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 


B 

3*5 

3-17 

5-17 

3 

5 

7.23*89 

599479 

7 

43691 

3 

71.122921 

127 

31 

31-127 

5.7.3.3.13.73 

9.7.5.73.37.109 
5.7.27.19.37.109 

5.7.27.19.13 

9.7.73.13.109.37 

3.3.7.5.73 

9.5.7.37.109 

9.7.13.37.109 

9.7.13.73 
3.5.7.13.73.19 

3-3.3. 5-7 
9.7.19.13 
9.5.13.19.73 

3.3.5. 7 
3-3.7. 13 

5.9.13.73 
3* 7-13- 19 

3.3.5.73 
3.3.73 

3.5.7.19 

7-13 

9*13 

3.3.7 
9-73 

73 

3.7*19 

5*3*3 

13 

3*3 

7 

5 

3 


e(B) 

4 

8 

8 

2 

4 
33 

33 
3 

34 
2 

35 
7 

5 

35 

36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
12 
12 
36 
36 
36 
18 
36 
12 
12 

6 

18 
9 

18 
12 
12 
6 

3 

4 
2 


D  (lower  bound) 


m 


3 

5 

5 
0 
1 

13 
18 

2 

14 
0 

22 

6 
4 

11 

17 

26 

25 

17 

27 

13 

19 

20 

14 

19 
9 

14 

20 

7 
9 

15 
11 
11 

8 

10 

5 

6 
5 
8 
5 
8 
4 
3 
2 
2 
1 
0 


8 

8 

8 

16 

16 

3 

3 

11 

4 
17 

3 

5 
7 
5 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 
3 

3 
6 
6 
6 

4 
4 
6 
6 
9 
9 
6 
6 
8 
6 
9 

12 

12 

12 

18 

18 


^(actual) 

8 

8 

8 

16 

16 

3 

3 

11 

4 
17 

4 

5 
7 
5 
3 
3 
3 
3 

3 

4 
4 
4 
4 
4 
4 
4 

4 
6 
6 
6 

5 
8 
8 
8 
9 
9 

6 
8 
8 
8 
9 

12 

12 

12 

18 

18 
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From  the  above  tabulated  results,  it  is  seen  that  the  lower 
bound  developed  in  this  section  is  very  tight.  This  lower  bound  also 
su86ests  a  criteria  for  one  to  choose  proper  B's  in  constructing  any  cyclic 
AN-code.  Suppose  that  we  want  to  construct  a  cyclic  AN-code  with  a  speci-  . 
fied  error  correcting  capability.  First  we  let  the  lower  bound  on  the 
code  distance  to  be  D^,  if  the  code  length  n  is  odd  and  the  code  is  so 
constructed  that  e(B)"n,  (other  cases  can  be  considered  in  a  similar  man¬ 
ner),  then  by  theorem  3.5.2,  we  have 


Di*  r— i+2-r^i  « 

Dj-1  =  (n-2)  /ij-u  v 

Dx-1  <  (n-2)/ix 
ix  <  (n-2)/(D1-l) 

4  i  ro^/CDj-iyi 

log2(B/3)-u’  1  f (n-2) / (Dj-iy) 


where  0<u<l 


where  0<u'<l 


the  smallest  integer  that  is  greater  than  or  equal  to  log2(B/3)  is  thus. 


hence. 


log2(B/3)  <  f(n-2)/(D1-l)l  +1 


B/3  <  2  f(n-2)/(Di-l)l  +1 


or, 

B  <3.2  r(n-2)/(D1-l)l  +1  A 

From  the  above  expression,  we  see  that  for  a  fixed  code  length  n,  the 
number  of  code  words  should  be  smaller  than  B^  so  that  the  code  constructed 
guarantees  to  possess  a  minimum  distance  at  least  as  great  as  D^.  Let  us 
consider  the  following  example: 


Suppose  that  we  want  to  construct  a  single  error  correcting 


cyclic  AN-code  of  length  15.  The  code  must  have  distance  at  least  3,  and 

can  be  calculated  as 

B  -  3-2  r(n-2)/<D1-l)T  +1 

.  3.2  f (15-2)/ (3-1)1  +1 

-  3‘27  -  384 

We  choose  B  to  be  smaller  than  384,  (e.g.,  B«31*7-217)  and  e(B)-15,  the 
code  is  of  minimum  distance  at  least  3. 

By  using  the  lower  bound  on  the  minimum  distance  derived  in 
theorem  3.5.2,  we  can  synthesize  any  cyclic  AN-code.  for  the  correction  of 
a  specified  number  of  errors.  The  following  codes  are  obtained: 


n  = 

length  of  the  code 

B  = 

number  of  code  words 

k  = 

number  of  information  bits 

D  .  - 
ml 

lower  bound 

of  the  minimum  distance 

D  = 
ma 

actual  minimum  distance 

r  = 

rate  of  the 

code 

n 

B 

k 

D  i 

D 

r 

ml 

ma 

42 

5419 

13 

6 

6 

0.31 

42 

16257 

14 

6 

6 

0.33 

44 

10565 

14 

6 

6 

0.32 

48 

24929 

15 

6 

6 

0.314 

50 

8283 

14 

8 

10 

0.28 

50 

44561 

16 

8 

8 

0.32 

52 

85489 

17 

8 

8 

0.33 

52 

253241 

18 

5 

6 

0.346 

52 

1266205 

21 

6 

6 

0.404 

54 

784899 

20 

6 

8 

0.37 

54 

87211 

17 

6 

6 

0.315 

54 

261633 

18 

8 

8 

0.33 

58 

3033169 

22 

5 

6 

0.38 

"fvrv 


429325 

3520465 

162565 

402905 

812825 

270805 

6700417 

1397419 

4192257 

411849 

62571 

45761 

652805 

130561 

3505429 

25080101 

12083 

36249 

398739 

3705353 

11116059 

16773121 

1774001 

9335617 

104323 

25781083 

77343249 


lilil — ^pper  Bound  on  the  Plat- an™ 


0.315 

0.365 

0.30 

0.317 

0.33 

0.317 

0.36 

0.318 

0.33 

0.288 

0.244 

0.244 

0.295 

0.25 

0.325 

0.368 

0.20 

0.228 

0.271 

0.315 

0.342 

0.33 

0.292 

0.33 

0.237 

0.338 

0.365 


To  find  an  upper  bound  on  the  distance  of  a  cyclic  AN-code, 
we  first  consider  the  congruence: 


2^  = 


mod  B 


(3.5.9) 
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I 

i 

I 

I 

i 

I 

I 

I 

I 

I 

I 

I 


where  0  £  i  <  n-1  is  given,  and  x  is  any  integer.  We  want  to  solve  for  the 
integer  N. 

Since  B  is  an  odd  integer,  2*  is  relatively  prime  with  respect  to  B 
for  all  is.  By  Euclidean  lemma,  we  have: 

s-21  +  t*B  =  1  for  some  integers  s  and  t  (3.5.10) 


Hence,  for  any  integer  x,  the  above  congreuence  has  a  unique  residue  class 
solution  for  N  modulo  B.  If  the  value  of  N  is  restricted  to  be:  OS  N<  B-l,  then 
the  solution  of  (3.5.9)  is  the  actual  value  of  N. 

For  a  code  word  AN,  0<  N<  B-l,  the  residue  2*N  mod  B  indicates  the 

nonzero  terms  in  its  NAF.  If  we  take  a  fixed  integer  i,  and  set  x  to  be  an 

integer  in  Mg,  then  the  solution  of  N  in  the  above  congruence  says  that  the  code 

word  AN  has  a  nonzero  term  in  its  NAF  at  the  position  2n  Also,  for  a  fixed  i, 

the  solution  N  of  the  congruence  (3.5.9)  are  all  distinct  for  different  x's  in  Mg; 
this  is  because 


2*N^  s  x  mod  B 
2^2  =  y  mod  B 

would  imply  that  x  s  y  mod  B  and  0<  x,y<  B-l  then  implies  x  «  y. 

These  set  of  N’s  obtained  by  choosing  a  fixed  i  and  setting  x  to  be  all 
the  integers  in  Mg  are  simply  the  codewords  that  have  nonzero  term  in  their  NAF's 
at  the  position  2  .  As  mentioned  above,  the  congruence  (3.5.9)  always  has 

unique  solutions.  That  means  we  are  always  able  to  find  the  code  words  with 
nonzero  terms  at  the  position  2n~i,  for  i  =  0,1,2,...,  n-1. 
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Fih-1  1 

There  are  ~T~  integers  in  the  region 
can  solve  an  N  with  respect  to  a  fixed  i.  In  other  words,  in  the  NAF  of  all  the 
code  words,  there  are  exact  code  words  with  nonzero  terms  at  the  position 

2n  This  argument  holds  true  for  all  i,  0^  i^  n-1,  hence  there  are  totally 

r&t-i  i 

n-  — —  nonzero,  terms  which  are  distributed  among  the  NAF's  of  all  the  nonzero 

codewords . 

The  number  of  the  nonzero  codewords  is  B-l,  hence  if  the  minimum 

distance  of  the  code  is  D  .  ,  we  have  the  following  theorem: 

mm 

Theorem  3.5.2:  The  minimum  distance,  D^n,  of  a  cyclic  AN-code  of  length 
n  is  at  most  as  large  as: 

Dmin  *  (B_1)  5  l^'1  *  n  (3.5.11) 

where  B  is  the  number  of  code  words. 

It  is  quite  interesting  to  note  that  for  large  B's,  this  bound  is 
approximately  equal  to  that  means  the  distance  of  the  code  can  not  exceed 
one  third  of  the  length  of  the  code! 

If  we  fix  the  length  n  of  the  code,  then  for  increasing  values  of  B, 
the  upper  bound  decreases.  This  is  because  the  ret.-?  of  the  code  increases  as  the 
number  of  code  words  increases,  while  the  error-c.  .ecting  capability  of  the  code 
decreases  as  the  rate  gets  higher. 


,  for  each  integer  in  we 


w*  i LV t  rf 
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4.  IMPLEMENTATION  OF  ARITHMETIC  CODES  BY  USING  MAJORITY  LOGIC 

4.1  Introduction 

One  of  the  major  unsolved  problems  in  the  theory  of  arithmetic  codes 
is  the  inplementation  of  decoding  circuitry.  Because  the  decoders  must  be  used 
frequently  in  systems  using  coding  redundancy,  the  practicality  of  coding  re¬ 
dundancy  depends  on  how  complex  the  decoders  are.  Also,  unlike  the  usual  situa¬ 
tion  in  a  data  transmission  system,  the  delay  introduced  in  the  decoding  circuit 
is  always  important  to  a  digital  computing  system.  Consequently,  a  practical 
decoder  should  have  simple  logic  circuitry  and  the  decoding  procedure  should  be 
fast. 

The  decoding  of  arithmetic  codes  for  multiple  error  correction  has 
been  shown  to  be  difficult.  Since  the  operation  of  the  codes  is  an  integer  or 
modulo  arithmetic,  the  only  information  that  decorder  receives  is  one  equation 
(  or  one  congruent  relation)  in  several  integral  unknowns.  In  terms  of  number 
theory,  the  decoding  equation  is  usually  a  Diophantine  equation,  which  is  extremely 
difficult  to  be  solved.  The  only  known  decoding  scheme  for  multiple  error  correc¬ 
tion  arithmetic  codes  is  the  "permutation  of  residue"  technique,  obtained  by 
Laste  and  Tsao-Wu  [30]  and  Hong  [29],  In  the  following,  we  shall  briefly  discuss 
the  permutation  of  residue  technique,  to  see  how  difficult  this  approach  is. 

The  first  step  of  the  "permutation  residue"  decoding  scheme  is  to 
calculate  the  decoding  index  k,  which  is  defined  as  the  smallest  positive  integer 
such  that  k  A  >  m®x  |e  .  I:  where  E  is  the  error  and  E  .  is  some  cyclic  shift  of 
E.  The  syndrome  S  s  E  mod  A  is  then  computed;  calculate  s  2  S  mod  A  for  each 
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left  shift  of  S  mod' A.  For  each  i,  the  arithmetic  weight  of  is  computed,  to 

see  if  W(S^)  is  within  the  error  correcting  capability  of  the  code.  If  none  of 

the  S^'s  satisfy  W(S^)  <  t,  where  t  is  the  error-correcting  capability  of  the 

code,  the  =  A  -  =  21(A-S)  mod  A  are  calculated  one  by  one;  for  each  i  the 

arithmetic  weight  W(Sp  is  counted  to  see  if  W(S|)  ^  t.  The  same  procedure 

goes  through  all  +  jA  and  S|  +  jA  for  all  0  <  j  <  k-1,  0<  i  <  n,  where  n  is 

the  length  of  the  code,  until  for  some  i^,  j^,  the  inequalities  W(S^  +  j^A)  :S  t 

d 

or  W(S|  +  j^A)  <  t  are  satisfied,  we  obtain  the  i^-th  cyclic  shift  of  the  error 
d 

pattern  E. 

It  is  easy  to  see  that  the  essence  of  this  approach  is  a  table  look-up 
attempt.  Although  the  amount  of  search  is  reduced  from  that  required  by  a  brute 
force  approach,  the  "permutation  of  residue"  scheme  still  involves  a  great  deal 
of  searching.  Also,  a  great  deal  of  complicated  calculations  must  be  done  by  the 
decoder,  and  the  decoding  delay  would  be  considerable  due  to  the  large  number  of 
shifting  cycles. 

If  the  "permutation  of  residue"  decoding  scheme  was  implemented  in  a 
digital  computing  system,  it  would  be  natural  for  one  to  ask  "can  this  decoder 
be  used  to  provide  as  good  or  better  use  of  redundant  parts  than  do  other 
techniques?"  Therefore,  how  to  find  decoding  algorithms  that  can  be  easily 
implemented  is  a  problem  of  both  theoretical  importance  and  practical  signifi¬ 
cance. 

Since  the  majority-logic  decoding  (MLD)  scheme  is  by  far  the  easiest 
one  that  can  be  applied  to  many  data  communication  channels,  it  would  certainly 
be  very  interesting  if  the  decoding  of  arithmetic  codes  in  a  digital  computing 
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system  can  also  be  done  by  using  majority  logic.  This  is  not  clear,  however, 
due  to  the  propagation  of  arithmetic  errors  as  a  result  of  carry  or  borrow 
failure. 

In  this  chapter,  we  are  able  to  divise  a  majority-logic  decoding 
scheme  for  some  arithmetic  AN-codes.  We  will  show  that  the  cyclic  AN-codes 
constructed  in  the  previous  chapter  are  all  majority-logic  decodable.  Further¬ 
more,  this  majority-logic  decoding  scheme  can  be  generalized  to  a  multi-step 
majority-logic  decoding  algorithm,  the  number  of  steps  required  for  decoding 
is  related  to  the  number  of  prime  factors  contained  in  the  length  of  the  code. 
With  this  majority-logic  decoding  scheme,  the  implementation  of  the  cyclic  AN- 
code  in  a  real  digital  computing  system  would  be  much  simpler. 

4.2  Errors  in  Modular  Arithmetic 

The  arithmetic  in  most  digital  computers  is  arithmetic  modulo  m  where 
either  ra  =  2n-l  (one’s  complement  arithmetic)  or  m  =  2n  (two’s  complement  arith¬ 
metic).  In  an  arithmetic  modulo  m,  the  integers  under  consideration  are  the  set 
of  integers  [0,  1,  2,  ...  ,  m-l},  and  the  additive  operation  is  defined  as: 

+  \  -  Km  (VV  (4.2.1) 

where  0  <  1^,  I  <  m-l,  and  Rm(x)  denotes  the  residue  of  x  modulo  m.  Similarly 
the  multiplicative  operation  is  defined  as, 

ri  •  J2  ■  \  <YV  (4-2-2) 

The  set  of  integers  {0,  1,  2,  ...,  m-l)  under  the  above  two  operations  forms  a 
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ring,  called  the  ring  of  integers  modulo  m.  For  any  integer  I  f  0  in  the  ring, 
its  negative  (or  the  additive  inverse)  is  the  integer  m-I. 

In  an  arithmetic  unit  which  performs  arithmetic  modulo  m,  we  shall 
always  assume  that  the  result  of  any  operation  is  represented  by  some  integer  R 
such  that  0  <  R  <  m.  If  the  arithmetic  is  modulo  2n-l,  we  have  an  n-stage  one's 
complement  adder.  A  cyclic  AN-code  of  length  n,  where  AB  =  2n-l,  can  be  used 
for  error-detection  or  correction  in  this  adder.  When  two  code  words  AN^  and 
AN^  are  added,  the  sum  has  the  form: 

Sum  of  (AN-.+AN-)  =  R  (AN..4AN„)  where  m  =  2n-l 

1  L  m  L  i 

-  A'R^  (Nj+N2)  (4.2.3) 

Since  0<  R^(Nj+N2)  <  B,  it  follows  that  the  modulo  sum  of  any  two  code  words  is 
another  code  word  in  the  cyclic  AN-code.  Hence,  a  cyclic  AN-code  is  linear  under 
modulo  arithmetic. 

If  an  error  has  occurred  during  the  additive  operation,  we  first  note 
that  the  absolute  value  of  the  error  is  less  than  2n-l,  for  the  error  is  caused 
by  adding  logic  failure  and  there  are  only  n  adding  units  in  the  adder.  An  error 
may  change  the  correct  sum  (which  is  a  code  word)  into  an  integer  which  is  greater 
than  2n-l,  in  that  case  the  actual  result  from  the  output  of  the  adder  is  the 
modulo  sum  of  the  correct  code  word  and  the  error  pattern: 

R  =  actual  sum  =  R  [A*R  (N.+N„)  +  E  ]  (4.2.4) 

m  Biz 

The  syndrome  S  associated  with  the  possibly  erroneous  result  R  is  de¬ 
fined  to  be  the  residue  of  R  modulo  A,  i.e.  S  *»  R^(R).  Since  2n-l  =  m  ■  AB, 
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we  have; 


s  =  Ka(R)  »  RA[RA3(A'RB(Ni+N2))  +  E  ] 

=  Ra  (E)  (4.2.5) 

Therefore,  the  syndrome  S  is  uniquely  determined  by  the  residue  of  E 
modulo  A.  There  are  B  distinct  error  patterns,  namely,  E,  E  +  A,  E  +  2A,  ...  , 

E  +  (B-l)  A  which  have  the  same  syndrome  S,  they  form  a  coset  expansion  of  the 
principal  ideal  (AN|0S  N<  B-l}  in  the  ring  of  integers  modulo  2n-l.  The  coset 
leader  is  the  pattern  that  has  a  minimum  arithmetic  weight;  it  is  the  error 
pattern  that  is  most  likely  to  occur.  The  d  v.  der  then  attempts  to  identify 
this  most  likely  occurred  error  pattern  from  the  syndrome  by  some  intellegent 
scheme. 

4.3  Majority  Logic  Decoding  Scheme  for  Arithmetic  Codes 

Let  us  consider  &  cyclic  AN-code  with  a  composite  block  length 

nl  n 

n  =  n^'n^.  Since  njn  implies  (2  -1)|(2  -1),  we  may  choose  the  generator 

A  as: 
o 

A0  =  (2n-l)/(2  X-l) 

n,(n  -1)  n.(n0-2)  n. 

=2  1  +  2  1  Z  +  ...  +  2  +  1  (4.3.1) 

In  this  A  N-code,  all  codewords  are  of  the  form  A  N  with  N  =  0,  1,  2,  ...  , 
o’  o  ’ 

nl 

2  -2.  If  we  divide  the  n  bits  of  a  codeword  into  n^  blocks,  then  each  block 

nl 

is  of  length  n^.  Since  N  is  less  than  2  -1,  we  have  enough  positions  to 
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express  N  in  its  binary  form  at  each  block.  Therefore,  the  binary  form  of  any 
codeword  AqN  consists  of  n2  replicas  of  the  binary  form  of  N.  It  has  been 
indicated  by  Erosh  [23]  that  the  minimum  distance  of  this  code  is  equal  to  n^, 
hence  this  code  can  correct  all  arithmetic  errors  of  weight  up  to  r(n2“l)/2l  . 

Suppose  that  two  coded  numbers  AqN^,  AqN2  are  added  in  an  adder  which 
performs  modulo  2n-l  arithmetic,  and  an  error  E  occurs,  then  the  incorrect  sum 
R  from  the  output  of  the  adder  is: 

R  5  Aq(N1+N2)  +  E  mod  2n-l  (4.3.2) 

Since  R  is  an  integer  in  the  ring  of  integers  modulo  2n-l, 

0<  R  <  2n-l;  we  have: 

R  s  A  N„  +  E  mod  2n-l  (4.3.3) 

o  3 

where 

K3  =  Rg  «,«2> 

We  assume  that  E  is  a  correctable  error  pattern,  then  W(E)  <  r(n2’l)/2l  . 

To  decode,  the  decoder  wants  to  recover  the  correct  sum  A  N_  from  the  res  .  R. 

7  o  3 

We  now  prove  the  following  lemma  and  theorems  which  yield  the  majority  decoding 
scheme: 


Lemma  4.3.1:  The  actual  sum  A  N„  is  a  zero  word  if  and  only  if  the  modular 

o  3 

weight  of  the  result  R  is  less  than  or  equal  to  [(n2~l)/2]. 

Proof:  Suppose  that  A^^  =  0,  the  (4.3.3)  can  be  written  as: 

R  =  E  mod  2n-l 


(4.3.4) 
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which  implies  that  R  =  E  for  E  >  0  or  R  =  (2n-l)  +  E  for  E  <  0.  Hence,  Wm(R) 

=  Wm(E)<[<l.2-1)/2]. 

On  the  other  hand,  if  Wffi(R)  <  [(n2-l)/2],  we  have: 

W(A  N  )  =  W  (R-E) 
m  o  3  m 

<  W  (R)  +  W  (E)  (by  triangle  inequality) 

mm  " 

<  C(n2-D2]  +  [(n2-l)/2] 

<  n  -1  (4.3.5) 

which  implies  that  AQNg  =  0  since  all  nonzero  code  words  are  of  weight  at 
least  n  . 

2  Q.E.D. 

Lemma  4.3.1:  determines  the  actual  sum  to  be  zero  or  not.  For  Wm(R)  > 
[(n2-l)/2],  we  have  the  following  arguments. 

We  first  consider  how  the  binary  bits  of  e  correct  sum  A^N^  are 
affected  by  a  single  arithmetic  error.  A  single  arithmetic  error  may  cause 
several  bits  to  be  in  error  because  of  the  influecne  of  carry  or  borrow  propa¬ 
gation.  Sometimes  an  error  may  cause  a  carrry  process  which  changes  the  code 
word  into  an  integer  greater  than  2n-l.  In  a  modulo  2n-l  adder  (i.e.,  one's 
complement  arithmetic),  there  are  only  n  positions  for  the  binary  representation 
of  the  result  R.  Hence,  the  carry  propagation  may  go  beyond  the  bit  with  position 
2n  *  and  then  change  cyclically  the  bits  with  positions  2^,  2^,  2^  ...  and  so  on. 
However,  any  carry  propagation  will  stop  whenever  a  zero  bit  is  reached. 

Similarly,  any  borrow  propagation  will  stop  whenever  a  bit  1  is  reached.  Since 
the  generator  A_^  of  the  code  is  of  the  form  as  shown  in  (4.3.1),  any  nonzero 
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code  word  can  have  at  most  n^-1  consecutive  0’s  or  1*8  in  its. binary  form  (we 
note  that  the  bits  with  positions  2n  ^  and  2^  are  considered  to  be  consecutive). 
Starting  at  the  erroneous  bit,  any  carry  or  borrow  propagation  caused  by  a  single 
arithmetic  error  cannot  propagate  more  than  n^  consecutive  positions.  In  other 
words,  the  carry  or  borrow  propagation  caused  by  a  single  arithmetic  error  can 
change  at  most  n^  consecutive  bits  of  a  nonzero  code  word  in  the  A^N-code.  Thus, 
we  have  proved  the  following  theorem: 

Theorem  4.3.1:  A  single  arithmetic  error  at  any  bit  of  a  nonzero  code  word 
in  the  A^N-code  can  change  at  most  n^  consecutive  bits  of  this  code  word. 

Next,  we  consider  a  multiple,  but  correctable  error  pattern  E.  Ex¬ 
pressing  the  error  E  in  its  NAF  as: 

~*1  —*2  *t 

E  =  ex  2  +  e2  2  +  . . .  +  et  2 

where  e^  *  +  1;  t<  [(n2-l)/2]  (4.3.6) 

We  partition  the  n  binary  bits  of  R  into  n^  disjoint  sets,  each  set 

.  n.+k  2n.+k  (n.-ljn^+k 

contains  n^  bits  with  positions  2,2  ,2  ,  ...  2  for  a  fixed  k, 

where  OS  k$  n^-1.  For  any  fixed  k,  these  n2  bits  are  ^  +  1  positions  apart 

from  one  another,  therefore,  by  theorem  4.3.1,  no  two  bits  of  the  same  set  can 

be  changed  by  a  single  arithmetic  error.  Since  the  arithmetic  weight  of  E  is 

k  nj+k  2nj+k  (n2-l)n1+k 

t,  at  most  t  of  the  n?  bits  with  positions  2,2  ,2  ,  ...  2 

can  be  altered  by  E,  the  rest  remain  unchanged  as  in  the  binary  form  of  the 

correct  sum  A  N. .  But  E  is  a  correctable  error,  t  is  less  than  or  equal  to 
o  3 

[(n2-l)/2],  thus  we  have  proved  the  following  theorem: 


Theorem  4.3.2:  If  the  error  is  of  weight  less  than  or  equal  to  [(n2"l)/2]> 
then  in  the  binary  form  of  R,  the  majority  of  the  bits  with  positions: 

.  n.+k  2n.+k  (n.-l)n  +k 

2K'2i2i  2  1  1 

where 

0  <  k  <  nL-l 

remain  the  same  as  in  the  binary  form  of  the  correct  code  word. 

From  the  lemma  and  theorems,  we  now  summarize  the  decoding  algorithm 

nl 

for  the  AQN-code  generated  by  Aq  =  (2n-I)/(2  -1)  as  follows: 

(1)  If  Wm(R)  <  C(n2-l)/2],  decode  AqN3  =  0. 

(2)  If  W  (R)  >  [(n_-l)/?.],  work  on  the  n1  sets  of  bits  in  the  binary 

m  "  .  n.+k  zn.+k  (n_-l)n..+k 

form  of  R  with  positions  2  ,  2  ,2  ,  ...,  2  for 

k  =  0,  1,  2,  ...,  n^-1.  For  each  set  take  the  majority  value  of 

the  bits,  then  form  a  block  of  n^  bits  with  the  majority  value 

obtained . 

It  is  straightforward  to  see  that  the  above  majority  logic  decoding 

n 

scheme  is  also  applicable  to  the  A'N-code  generated  by  A*  =  (2  - 1) / (2  -1).  The 

v  O 

decoding  procedure  is  exactly  the  same  except  we  divide  the  n  bits  of  R  into  n^ 
blocks,  with  each  block  of  length 

Since  there  is  no  restriction  on  the  parameters  n^  and  n^,  the  majority 

logic  decoding  scheme  can  be  applied  quite  generally  to  any  AqN  or  A^N-codes  of 

composite  length.  When  n^  and  n2  are  primes,  it  should  be  mentioned  that  the 

A^N-code  is  the  p^p2~code  generated  by  A^  and  the  A^N-code  is  the  p^p2~code 

generated  by  A  ,  as  discussed  in  Chapter  3. 

6 
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4.4  2-Step  Majority-Logic  Decoding  Scheme 

It  has  been  shown  in  the  previous  section  that  the  cyclic  AoN-code  is 

majority- logic  decodable.  To  decode,  we  only  need  one  level  of  majority  logic. 

It  is  natural  to  that  by  using  more  than  one  level  of  majority  logic  other 

cyclic  AN-codes  might  also  be  decoded  easily.  In  this  section,  we  will  derive 

a  2-step  majority- logic  decoding  scheme  which  can  be  applied  to  the  class  of 

multiple  error-correcting  AN-code  described  in  Chapter  3. 

We  recall  that  the  codes  described  in  Section  3-.3are  generated  by  the 

nl  n2 

cyclotomic  generator  A  =  (2  -l)/(2  -1)(2  -1),  where  (n^n^)  =  1  and  n  =  n^n,,. 

If  n^  and  n^  are  primes,  then  the  codes  are  the  p^p^-codes  generated  by  A^  in 

Section  3.4.  If  we  impose  one  additional  condition  on  the  parameters  n^  and  n^ 

so  that  n^  ^  2n^-l,  then  the  majority- logic  decoding  can  be  done  as  follows: 

From  theorem  3.3.1,  we  know  that  the  minimum  distance  of  the  cyclic 
n  nl  n2 

AN-code  generated  by  A  =  (2  -l)/(2  - 1)  (2  -1)  is  n^  provided  that  (n  .n^  =  1 

and  n^  <  n^.  Therefore,  a  correctable  error  pattern  E  is  of  weight  W(E)  < 

C (tt1“l)/2] . 

Suppose  that  two  coded  numbers  AN^  and  AN^  are  added  in  an  adder  which 
performs  modulo  2n-l  arithmetic  and  a.  correctable  error  E  occurs,  then  the  in¬ 
correct  sum  R  from  the  output  of  the  adder  is: 

R  s  AN3  +  E  mod  2n-l, 

where 

n3  -  R^-rty 

W(E)  <  [(n^D/2] 


(4.4.1) 
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I 

I 

I 

I 

I 

I 

I 

I 


For  simplicity,  we  denote  2-1  by  m  in  the  following  discussion. 

n2 

Multiplying  both  sides  of  (4.4.1)  by  (2  -1),  we  have: 

R(2  -1)  s  A (2  -1)  •  N3  +  E(2  -1) 


(2n-l) 


nl  n2 
(2  -1)(2  -1) 

n„ 


n„  n„ 

•  (2  -1)  N,  +  E(2  -1) 


s  A  N,  +  E(2  2-l)  mod  2n-l 
o  3 


(4.4.2) 


where  Aq  =  (2  -l)/(2  -1).  Denoting  the  residue  of  any  integer  I  modulo  2  -1 

by  R  (I),  then  R  [A  N„1  is  a  code  word  in  the  A  N-code  discussed  in  previous 
m  m  o  3J  o  „  „ 

°2  nln2 

sections.  Expressing  E(2  -1)  in  NAF  and  substituting  1  for  2  ,  we  can 

n2 

reduce  E(2  -1)  to  an  integer  E  with  the  property  that  the  absolute  value  of 

n  K 

E  ,  i.e.  |E  I,  is  less  than  2  -1.  We  note  that  E,.  is  not  necessarily  the 

TV'  l 

n  n 

same  as  R  (E(2  -1))  since  R  (E (2  -I))  is  always  prositive  while  E  may  be 
m  m  r 


negative. 


Equation  (4.4.2)  can  now  be  rewritten  as: 


R  [R(2  -1)]  =  R  [R  (A  N,)  +  E  ] 

nr  v  rav  o  3  rJ 


(4.4.3) 


The  integer  R  [R(2  -1)]  can  be  considered  as  a  corrupted  code  word 


m 


of  the  code  generated  by  Aq.  Since  the  cyclic  AQN-coae  is  of  distance  n^,  and: 


) 

•  *■ 


W(E)  <:  [(n^D/2] 


n2  2nrx 


I 

I 


mm 
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W(Er)  <S  W(E(2  2-l)) 

<  2W(E)  (by  triangle  inequality) 

<  nrl<  [(n2-l)/2] 

the  error  pattern  E  is  within  the  error-correcting  capability  of  the  cyclic 

n2 

A^N-code.  Therefore  by  div  ..  i-*  the  n  bits  of  Rm[R(2  -1)]  into  n2  blocks  the 

majority- logic  decoding  scheme  for  the  AQN-code  will  correctly  yield  the  code 

word  R  (A  N  ), 
m  o  3 

To  obtain  the  actual  error  pattern  E,  we  need  the  second  step.  Denoting 

n? 

{Rm[R(2  -1)]  -  R^^N^)}  by  E',  it  can  be  seen  from  equations  (4.4.2)  and  (4.4.3) 

that: 

n 

E*  s  E(2  2-l)  mod  2n-l  (4.4.4) 

hence, 

n,  n  n9 

E* / (2  -1)  5  E  mod  (2n- 1) / (2  -1)  (4.4.5) 

n  n? 

Recalling  that  (2-1)/ (2  -1)  =  A'q,  we  have: 

n9 

Rm[E'/(2  -1)]  =  RjR/A’  *  k)  +  E  ]  for  some  k  (4.4.6) 


1 


n2 

From  (4.4.6),  the  integer  Rm[E'/(2  -1)]  can  be  considered  as  a 

corrupted  code  word  of  the  cyclic  AN-code  generated  by  A^.  Since  the  cyclic 

A^N-code  is  of  minimum  distance  n^,  and  since  W(E)  ;£  £(n^-l)/2],  the  actual 

error  pattern  E  can  be  decoded  by  the  majority-logic  scheme  for  the  cyclic 

A'N-code. 

o 


jiia 


The  2-step  majority-logic  decoding  scheme  derived  above  can  now  be 


summarized  as  follows: 

(1)  If  W(R)  S  [(nj-l)/2],  decode  AN^  =  0. 

(2)  If  W(R)  >  [ (n^-l)/2},  get  E*  by  applying  the  majority- logic  decoding 

n2 

scheme  for  AoN-code  to  R^[R(2  -1)].  (first  step) 

(3)  Obtain  the  actual  error  pattern  E  by  applying  the  majority- logic 

n2 

decoding  scheme  for  A^N-code  to  Rm[E'/(2  -1)3 -  (second  step) 

It  should  be  remarked  that  whilq  the  condition  n  ^  2n^-l  *mPose<*  on 

nln2  ni  n? 

the  codes  generated  by  A  =  (2  -l)/(2  -1)(2  -1)  guarantees  the  corection  of 

errors  up  to  minimum  distance,  the  condition  is  not  necessary  for  the  2-step 

majority- logic  decoding  scheme  to  work.  In  general  the  majority- logic  decoding 

scheme  will  correct  up  to  error  patterns  of  arithmetic  weight  the  smaller  of 

C(nrl)/2]  and  \  [(n2-l)/2]. 

Let  us  consider  the  following  example: 

Example:  Suppose  that  the  block  length  of  the  code  n  =  5.9  =  45,  and  the 
45  9  5 

generator  A  =  (2  - 1) / (2  -1)(2  -1).  Then  the  code  is  of  minimum  distance  5, 

capable  of  correcting  any  double  arithmetic  errors.  We  have. 


A  -  -  =  231  +  226  +  222  +  ,21  •  l17  +  216  f  213  +  212  *  2U 

(2V-1) (2^-1) 

+  2+2  +  2"  +  2  ' 


6  /•  3  2 

' v  *  «  *  j  v.  '  .  r.  ** 


+  2+1 


=  (000000000000010000100011000110011100111011111) 


If  two  coded  numbers  37A  and  27A  are  added  in  a  45-stage  one’s  comple^ 

38  8 

ment  adder,  and  a  double  error  E  = *2  +2  occurs,  then  the  incorrect  sum  R  from 

the  output  of  the  adder  is,  R  =  A(37+27)  +  E  =  64A  +  E.  To  decode,  we  first 
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9  45 

calculate  the  residue  o £  R(2  -1)  modulo  2  -1  as: 

00010,  01010,  00010,  00010,  00010,  00110,  00001,  11001,  11110 
We  divide  the  above  binary  form  into  9  blocks  and  check  that. 


digits  with  positons  majority  value 

5k+4 

2  ,  0<S  kS  8  0,  0,  0,  0,  0,  0,  0,  1,  1  0 

5k+^ 

2  ,  0<  k  <  8  0,  1,  0,  0,  0,  0,  0,  1,  1  0 

5k+2 

2  ,  0<  8  0,  0,  0,  0,  0,  1,  0,  0,  1  0 

5k+l 

2  ,  0<  k  <  8  1,  1,  1,  1,  1,  1,  0,  0,  1  1 

25k+0,  0<  k^  8  0,  0,  0,  0,  0,  0,  1,  1,  0  0 


q 

Hence  the  majority  decision  on  Rm(R(2  -1))  yeidls  a  code  word: 

Rm(AoN)  »  00010,  00010,  00010,  00010,  00010,  00010,  00010,  00010,  00010 

The  error  at  this  state  is  E*  =  R  (R(29-l))  -  R  (AN),  which  has  the 

m  mo 

binary  form: 


00000,  01000,  00000,  00000,  00000,  00100,  000(-l)l,  110(-1)1,  11100 

9  45 

The  actual  error  pattern  F,  is  congruent  to  E'/(2  -1)  modulo  2  -1, 


E(29-l)  s  E'  =  238  +  21?  -  28  -  22 


mod  2  - 1 


E  s  229  +  220  +  211  +  28  +  22  mod(245-l)/(29-l) 


which  has  the  binary  form: 


000000000,  000000100,  000000100,  000000100,  100000100 


(4.4.7) 
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Again,  the  majority-logic  scheme  on  (4.4.7)  yields  a  block  000000100. 

Repeating  this  block  five  times,  we  have: 

000000100,  000000100,  000000100,  000000100,  000000100  (4.4.8) 

The  binary  integer  (4.4.8)  is  a  code  word  generated  by  °»  (2^-1)  /  (2^-1) . 

Subtracting  (4.4.8)  from  (4.4.7),  we  get  the  actual  error  pattern  E: 

000000 (-1)00,  000000000,  000000000,  000000000,  100000000 

OQ  O 

Hence,  the  error  pattern  E  is  decoded  as  -2  +2  . 

4.5  L-Step  Majority-Logic  Decoding  Scheme  for  Arithmetic  Codes 

In  the  previous  section  we  have  sho<  n  that  the  cyclic  AN-code 
nln2  ■  nl  no 

generated  by  A  =■  (2  -l)/(2  -1)(2  -1)  is  2-step  majority-logic  decod- 

able  provided  that  n2  >_ 2n^-l,  If  n^  and  are  two  prime  integers,  then 
the  code  is  the  p^p^-code  generated  by  Ay.  The  reason  that  the  decoding 
needs  2-step  of  majority  logic  is  that  the  denominator  of  the  generator  A 
contains  two  factors  of  the  form  2a-l.  We  must  multiply  the  form  genera¬ 
tor  A  by  a  factor  2a-l  so  that  the  modified  generator  of  the  code  has  a 
form  the  same  as  Aq  or  A^,  then  the  majority- logic  decoding  scheme  can  be 
applied.  After  the  first  step  of  majority-logic  decision,  we  divide  the 
result  by  the  factor  2a-l,  the  generator  of  the  code  at  this  stage  is 

(2  — 1) / (2  -1),  which  has  the  same  form  as  A  or  A',  hence  the  second 

o  o 

step  of  majority  logic  is  applied. 

It  is  not  difficult  for  one  to  realize  that  the  number  of  factors  of 
the  form  2  -1  contained  in  the  denominator  of  the  generator  is  related  to  the 
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number  of  steps  required  for  majority-logic  decoding.  As  a  matter  of  fact,  the 

2-step  majority- logic  decoding  scheme  developed  in  the  previous  section  can  be 

generalized  in  a  straightforward  manner  to  the  following  theorem: 

Theorem  4.5.1:  Let  n^,  n^,  ...,  n^  be  L  pairwise  relatively  prime  integers 

w4 th  the  condition  that  n^+^  s  2n^  for  i  =  1,2,...,  L-l;  then  the  cyclic  AN-code 
L 

of  length  n  =  tt  n.  generated  by: 


L  n 

A  =  (2n-l)/  tt  (2  1-1) 


(4.5.1) 


can  be  majority  logically  decoded  in  exactly  L  stops. 

Proof:  The  proof  of  this  theorem  is  essentially  the  same  as  the  proof  of 
2-step  majority- logic  decoding  scheme.  We  start  with  the  decoding  equation: 


R  «  AN  +  E 


(4.5.2) 


L  ni 

Multiplying  both  sides  of  (4.5.2)  by  Tt  (2  -1),  we  have: 

i=2 

L  n.  L  n  L  n. 

R  ‘  tt  O  -1)  =  [A-  tt  (2  -1)]  N  +  E  •  tt  (2  M) 

i=2  i=2  i=2 


where 


L  n 

“  A  N  +  E  •  tt  (2  -1) 

i=2 


\  =  (2n-I)/(2  :-l) 


(4.5.3) 


Taking  the  residue  of  both  sides  of  (4.5.3)  modulo  2  -1.  we  have: 


L  n. 


RJR-  <2  -1)!  -  +  Enl 

1=2 


(4.5.4) 


«»*****»«! i&ft. 


L  n. 

Where  E  ,  is  obtained  by  substituting  1  for  2n  in  the  8i&  of  E*  tt,  (2  -1). 

T  i=2 

We  assume  that  the  error  pattern  E  is  of  weight  less  than  or  equal  to 

r  L“L  i 

(  tt  -l)/2  ,  where  Txl  denotes  the  integeral  part  of  x,  then  we  have: 
i=i  1 

W(E  , )  £  W[E-  n  (2^1)]  <S  2L_1  W(E)  <  21"1  r(LTT1n.-l)/21  (4.5.5) 

ri  i=2 '  i=l  1 


Therefore, 


ni+1  S  2n^  for  i=  1,  2 . L-l 


,L-l  L’1  ^  L 

2  *  rr  n.  <  tt  n. 

i=l  i=2  1 


(4.5.6) 


(4.5.7) 


4 


W(E  ,)<  ’{(  tt  n.)  -l)/2 
i=2  x 


(4.5.8) 


Since  the  integer  is  of  the  same  form  as  Aq  in  section  4.3,  it 

L 

generates  a  cyclic  A^N-code  of  minimum  distance  tt  n^.  Also,  R_m(A^N)  is  a 


L  n. 

code  word  in  the  A^N-code,  and  the  integer  Rm[R'  tt  (2  -1)]  can  be  considered 

as  a  corrupted  code  word  in  the  A^N-code,  the  error  E  ^  is  within  the  error 

correcting  capability  of  the  A^N-code.  Therefore,  the  majority- logic  decoding 

L  n- 

scheme  can  be  applied  to  R  [R*  tt  (2  1- 1)] ,  we  get  R  (A.N)  correctly. 

m  i=2  ra  1 

L  \ 

Defining  by  the  difference  of  R^[R*  tt  (2  -1)]  and  Rm(A^N)  we  have: 

L  n 

E  £  E-  tt  (2  “-l)  mod  2n-l 
i=2 


I- 


m 
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L  n. 


E  /(2  2-l)  s  E*  tt  <2  i-l)  mod  (2n-l)/(2  2-l)  (4.5.9) 

i=3 


R[E./(2  2-l)]  =  K  [R  (A  *t)  +  E  j  for  some  t  (4.5.10) 

mi  mm/  r  z 

n 

where  A„  =  (2n-l)/(2  -1)  and  E„  is  obtained  by  substituting  1  for  2n  in  the 

^  r/ 

L  n. 

NAF  of  E-  tt  (2  -1). 

i=3 

The  intege-  Rm[Ej/(2  ^-1)]  can  be  considered  as  i  corrupted  code  word 

^  L 

of  the  cyclic  AN-code  generated  by  A.,  which  is  of  minimum  distance  n. •  rr  n.  . 

1  i=3  1 

Since, 


L  n. 


W(E  )  <  W[E-  rr  (2  L-l)]  <  2L‘2  W(E)  <  ?L'2  f(  n  n  -lj/21 
i=3  i=l  1 


0L-2  ,L-2 

2  ‘  "  "i  =  nl  "  nl  *  tt  n 

i=l  i=2  1  i=3 


(4.5.11) 


r  h  i 

W(E  )  <  ((n.  *  tt  n  )-l)/2 


(4.5.12) 


the  error  Er2  is  a  correctable  pattern  with  respect  to  the  AN-code  generated  by 

n2 

A- .  Using  majority- logic  decoding  scheme  on  R  [E.1//2  -1)1  we  can  determine 

A  m  l 

n9 

Rm(A2-t)  correctly.  Defining  E2  by  the  difference  of  Rj.E^  (2  -1)]  and  Rm(A-t), 


we  have: 


L  n. 

E  s  E  •  tt  (2  -1)  mcd  2°-l 

i=3 


(4.5.13) 


ivw  TO  1  W*W<4!»  WKfBBBJH 
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which  has  the  same  form  as  the  congruence  (4.5.8).  The  same  majority  decoding 
method  can  be  applied  and  we  get  another  congruence  of  the  same  form.  So  in 
general,  at  the  h-th  stage  we  have: 

L  n. 


E,  1  s  E  •  n  (2  i-l)  mod  2n-l 
n  i=h 


n  L  n. 

E.  ,/( 2  -1)  =  E  •  it  (2  -1)  mod  A. 

h-i  i=n+l  “ 


(4.5.14) 


(4.5.15) 


where  -  (2n-l)/(2  M-l). 

Therefore,  at  the  L-th  stage, 


T 

E  s  E  •  (2  -1)  mod  2n-l 


E k_,/ (2  -1)  s  E  mod 

r  L'1  i 

We  see  that  W(E)  <  ((  tt  n.)-i)/2  and  the  cyclic  AN  code  generated  by 

i=l 

L-l 

A.  is  of  minimum  distance  tt  n. »  hence  the  majority  logic  decoding  shceme  on 
^  1  !■ 


i=l 


E^_1/(2  -1)  will  correctly  yield  the  actual  error  pattern  E. 


Q.E.D. 


In  the  proof  of  the  above  theorem,  we  have  assumed  that  the  weight  of 

r  L_1  i 

E  is  less  than  or  equal  to  ((  tt  n,.)-l)/2  ;  the  fact  that  the  L-step  majority- 

i=l  * 

logic  decoding  r-'n  correct  errors  up  to  this  weight  gives  the  lower  bound  on  the 


minimum  distance  ■  the  code,  hence  we  have  the  following: 


Corollary  4,5.1:  The  cyclic  AN-code  as  stated  in  the  theorem  4.5.1  is  of 

L-l 

distance  at  least  as  great  as  rr  n. . 

i=l  1 

The  result  of  the  above  theorem  is  a  straightforward  generalization  of 
the  2-step  majority- logic  decoding  scheme.  The  code  constructed  has  rather  high 
redundancy.  To  improve  the  efficiency  of  the  code,  we  want  to  reduce  some  un¬ 
necessary  factors  contained  in  A.  In  the  following,  we  shall  construct  another 
class  of  codes  which  is  also  L-step  majority-logic  decodable.  This  code  is 
generated  by  a  cyclotomic  generator,  it  has  better  efficiency  and  a  very  general 
structure.  Although  the  decoding  for  this  new  code  is  similar  to  that  of  the 
code  in  theorem  4.5.1,  we  shall  go  through  the  discussion  in  detail  because  some 
important  remarks  must  be  made. 

Definition  4.5.1:  Any  positive  integer  n  can  be  expressed  as  a  product  of 

Si 

powers  of  prime  factors,  i.e.  n  =  tt  P^  .  This  product  of  prime  powers  form 

i 

is  unique  for  any  n  and  is  called  the  canonical  form  of  n. 

Theorem  4.5.2:  Let  the  canonical  form  of  a  positive  integer  n  be 
L  Si 

n  =  tt  p.  .If  the  distinct  prime  factors  p.'s  satisfy  the  condition  that: 
i=l  1  1  ' 

p^  ^  2*"  *  (p^-1)  +  1  for  i  =  2,  3,  ...»  L;  then  the  cyclic  AN-code  of  block 
length  n  generated  by  Qn(2),  the  primitive  cyclotomic  factor,  can  be  majority- 
logic  decoded  in  exactly  L  steps. 


Proof:  First,  let  us  consider  the  minimum  distance  of  the  code.  Since 
n  n  n/p^ 

A  =  Q  (2),  the  integer  (2  -l)/(2  -1)  is  a  multiple  of  A,  hence  it  is  a  code 

n  n/pl 

word.  The  arithmetic  weight  of  (2  —  1 )/ (2  -1)  is  p^,  hence  the  minimum  dis¬ 

tance  of  the  code  is  at  most  p^.  If  E  is  a  correctable  error  pattern,  then 


W(E)  <  (p^-l)/2.  We  start  with  the  decoding  equation: 
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.R  =  AN  +  E 


(4.5.16) 


L-l  n/p 

Multiplying  both  sides  of  (4.5.16)  by  rr  (2  -1), 

i=l 


L-l  n/p  L-l  n/p  L-l  n/p 

R  •  rr  (2  -1)  =  A  •  tt  (2  -1)  .N  +  E-  tt  (2  -1)  (4.5.17) 


i  ,'  .  n/p 

Let  A1  =  A-k  =  (2n-l)/(2  *■-!;,  then  B  =  (2n-l)/A  ■  k*(2  -1)  and  the  code 

generated  by  A  is  of  minimum  distance  pt .  Since, 

1 


B  *  Qn(2)  =  (2n-l)  =  rr  Qd(2) 


(4.5.18) 


the  cyclotomic  factors  contained  in  B  are  all  the  Q  (2)'s  with  d  a  proper 

d  L  n/p^ 

divisor  of  a.  Clearly,  these  Q  (2)'s  are  also  contained  in  tt  (2  -1),  hence 

i=l 


L  "  n/pt  L-l  n/pi 

B  divides  tt  (2  -1);  which  implies  that  k  divides  tt  (2  -1). 

i=l  i=l 


Therefore,  (4.5.17)  can  be  rewritten  as: 

L-l  n/p.  L-l  n/p. 

R-  tt  (2  X-l)  =  (A-k)-M  +  E-  it  (2  x-l)  for  some  M 
i=l  i=l 


L-l  n/p. 

AM  +  E*  n  (2  1-1) 

L  i=l 


(4.5.19) 


Taking  the  residues  of  both  sides  of  (4.5.19)  modulo  2  -1,  we  have: 


L-l  n/p. 

RjR-  TT  (2  M)]  =  R  [R(A.M)  +  E  ,  ] 

m  J  m  m  i  rl 


(4.5.20) 


»s?y  iwrwsn/f  > 


where  Erl  is  obtained  by  substituting  1  for  2°  in  the  NAF  of  E-  n  . 

In  addition, 

pl  * 2L_1  <prl>  + 1 
WOO  <  (prl)/2 

L-l  n/p 

W<Erl)  <:  W[E*  tr  (2  *-1)]  S  2L_1  W(E) 

X*  X 

s  2*“  (pl-1)/2  (4.5.21) 

Hence,  Erl  is  a  correctable  error  pattern  with  respect  to  the  cyclic  AN-code 
,  L-l  n/p. 

generate  y  A^  Also,  Rm[R*  ^(2  -i)]  can  be  considered  as  a  corrupted  code 

word  of  the  A^N-code,  then  the  majority- logic  decoding  scheme  for  A,N-code  can 

be  applied  to  ^(2  l-l)]  and  the  code  word  Rm(A^M)  can  be  determined 

correctly.  Defining  *  by  the  difference  of  (2n/Pi-l)-|  and  R  (A  M) 

we  have:  «  "  * 


L-l  n/p. 

E1  5  2  '  tt  (2  1-l)  mocj  2n-l 

i=l 


(4.5.22) 


»  y,9n/pL-l  ,,  L-2  n/p  n/p. 

21/(2  -1)  s  E*  n  (2  -1)  mod  (2  -1) / (2  L_1- 

i=l 


n/PT-l 

Rn£El^2  =  Rro^Rm^A2,t:^  +  Er?3  for  some  t 


(4.5.23) 


(4.5.24) 


where  A,,  -  (2n-l)/ (2  L"1-l)  and  E^  is  obtained  by  substituting  1  for  2n 

the  NAF  of  E-  rr  (2  Pi-1).  since. 


1  — Iff"!-) »l»-W H1  >* ■' l»».f 


W(Er2)  <  W(E 
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L-2  n/p  L_„ 

it  (2  -1))  <  2  W(E) 

i=l 

<  2L"2 (p, - 1) / 2  <  (p.  ,-l)/2 

1  W  (4.5.25) 


2 


and  the  code  generated  by  A„  is  of  distance  Pj_  E  „  is  a  correctable  error. 

‘  "  ^PL-1 

Applying  majority-logic  decoding  scheme  to  R  [E./(2  -1)1,  we  can  determine 

n/pL_i 

the  code  word  R  (A2t)  correctly.  Defining  E2  by  the  difference  of  E^/ (2  -1) 

and  R  (A  *t),  we  have: 
m  2 

L'2  n/pi  n 

E  a  E-  tt  (2  -1)  mod  2  -1  (4.5.26) 

2  i=l 


which  has  the  same  form  as  the  congruence  (4.5.22).  The  same  majority  decoding 
method  can  be  applied  to  (4.5.26)  and  we  obtain  another  congruence  of  the  same 
form.  So  in  general,  at  the  h-th  stage  we  have: 

L-h+1  n/p 

E,  .  =  E*  TT  (2  -1)  mod  2  -1 

h-1  r=l 


Eh-i/<2 


n/p 


L-h+1 


-1)  5  E- 


L-h  n/p. 

n  (2  x-l) 
i=l 


mod 


*h 


wher  Ah  =  (2n-l)/(2  /PL'h+1_i). 

Therefore,  at  the  L-th  stage  we  have: 

n/pl  n 

Et  .  =  E  •  (2  -1)  mod  2-1 

u*  X 

n/p, 

El_1/(2  -1)  s  E  mod  Al 

n/P, 

Rm[E^  ^/ (2  -1)]  =  R^R^A^’u)  +  E]  i°r  some  u 


(4.5.27) 

(4.5.28) 


(4.5.29) 


where  =  (2  -l)/(2  -1).  Since  the  code 

tance  p,  and  W( E)  <  (p1-l)/2,  the  majority-logic  decoding  shceme  on 
n/Px 

R  [E.  ,/(2  -1)]  will  correctly  yield  the  actual  error  pattern  E. 

Q.E.D. 

As  mentioned  previously,  the  code  generated  by  Qn (2)  is  of  minimum 
distance  at  most  p^.  But  from  the  above  theorem,  it' has  been  shown  that  the 
code  can  correct  all  arithmetic  errors  of  weight  up  to  (p^-l)/2,  hence  we  have 
the  following: 

1*  s .  .  , 

Corollary  4.5.2:  If  the  code  length  n  =  rr  p^  X,  and  p^  s  21  (p^-1)  +  1 

i=l 

for  i  =  2,  3,  . . .  L,  s  £  0  for  all  i  =  1,  2,  ...,  L,  then  the  minimum  distance 
of  the  cyclic  AN-code  generated  by  Qn(2)  is  exactly  equal  to  p^. 

Again,  it  should  be  mentioned  that  while  the  conditions  p^  £  2*  1 
(p^-1)  +  1  for  i  =  2,  3,  . . .  L,  imposed  on  the  code  guarantee  the  correction 
of  errors  up  to  minimum  distance,  the  conditions  are  not  necessary  for  the 
L-step  majority-logic  decoding  scheme  to  work. 

The  majority-logic  decoding  scheme  developed  in  this  chapter  is  a  com¬ 
pletely  different  approach  as  compared  with  the  permutation  of  residues  technique, 
no  search  is  required  at  all.  Those  majoriy-logic  decodable  AN-codes  can  then 
be  implemented  easily.  It  is  our  hope  that  the  majority-logic  decoding  scheme 
presented  here  will  be  applicable  to  other  classes  of  AN-codes  and  will  lead 
to  other  practical  decoding  implementations. 


generated  by  is  of  minimum  dis- 
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5.  APPLICATIONS  OF  CYCLIC  AN-CODES  TO  COMPUTER  ARITHMETIC 


5.1  Concurrent  Checking  in  Computer  Arithmetic 

A  logic  fault  occurs  in  &  digital  computing  system  when  a  defective 

9 

component  or  external  interference  causes  a  deviation  of  a  logic  component  from 
its  prescribed  value.  In  general,  there  are  two  approaches  which  can  be  employed 
to  check  the  presence  of  errors.  In  periodic  checking,  ordinary  operation  is 
periodically  interrupted  and  a  checking  program  is  carried  out  which  detects  the 
presence  of  an  error  and  indicates  its  location.  In  concurrent  checking,  special 
logic  circuits  are  used  to  detect  the  presence  of  errors  in  computer  words  con¬ 
currently  with  ordinary  operation  of  the  computer.  Both  the  well-known  repli¬ 
cation  of  processors  with  outpuc  voting  and  the  arithmetic  encoding  of  the 
operands  followed  by  the  application  of  decoding  algorithm  to  the  results  are 
the  methods  of  concurrent  checking  in  digita1  computing  systems. 

In  the  triplication  of  processor  method,  the  majority  vote- taker  gives 
the  correct  result  only  when  two  or  three  of  the  identical  processors  operate 
normally.  If  more  than  one  processor  had  errors,  the  majority  voting  scheme 
would  not  work.  In  other  words,  the  application  of  the  triplication  of  processor 
is  based  on  the  assumption  that  all  the.  logic  faults  are  restricted  to  only  one 
of  the  three  identical  processors.  This  model  is  c  stionable  sometimes. 


especially  when  multiple  error  patterns  occur  in  a  complicated  processor. 

On  the.  other  hand,  when  the  arithmetic  coding  method  is  used,  the 


redundancy  is  incorporated  into  the  operands  themselves  which  are  being  pro¬ 
cessed.  A  multiple  error-correcting  code  is  necessary  for  some  processors  when 
single  error  protection  is  not  reliable  enough.  As  pointed  out  in  the  previous 
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chapters,  the  decoding  of  multiple  errors  is  much  more  difficult  than  the  decodin 
of  single  error,  the  reliability  of  the  complicated  decoder  circuitry  can  not  be 
overlooked.  The  assumption  that  decoders  are  error-free  then  is  not  reasonable. 

In  this  chapter,  a  new  concurrent  checking  technique,  called  the 
"separate  checking"  is  studied.  This  separate  checking  technique  possesses  the 
advantages  of  both  the  coding  approach  and  the  replication  of  processors.  It 
combines  the  ideas  of  coding  and  replication  of  processors  properly  so  that  the 
error  model  is  more  realistic  and  the  checking  procedure  is  much  simpler. 
Furthermore,  since  only  single  error-correcting  codes  are  used  in  this  techni¬ 
que,  the  decoding  is  simple  and  thus  reliable. 

The  application  of  cyclic  AN-codes  to  the  high  speed  multiplier  is 
also  investigated  in  this  chapter.  It  will  be  shown  that  the  majority-logic 
decodable  cyclic  AN-codes  constructed  in  the  previous  chapters  can  be  applied 
practically.  All  tnese  applications  render  the  evidence  that  the  arithmetic 
coding  is  a  promising  approach  to  the  study  of  digital  computing  systems  diag¬ 
nosis  and  reliability. 

5.2  Separate  Concurrent  Checking  Techniques 

Suppose  that  two  integers  N.  =  (a  .,  a  ,  ...,  a  ,  an)  and 

i  n-i  n-2  i  u 

N9  =  (bn^,  b^  2>  •  ••>  h^,  b^)  are  added  in  a  n-stage  one's  complement  binary 
adder.  If  an  error  E  occurs,  then  the  result  from  the  output  of  the  adder  is 

R  =  Rss[Res(N1+N2>  +  Ej 

where  Rest”)  denotes  the  residue  of  r.  modulo  2n-l. 


(5.2.1) 
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Since  the  operands  and  ^  are  not  coded,  it  is  not  possible 

for  one  to  correct  the  error  from  the  erroneous  sum  R.  Instead  of  putting  ■ 

redundancy  into  and  N^,  we  use  several  separate  checkers  to  execute  the 

fault-tolerant  operation  as  follows: 

Let  us  call  the  n-stage  one’s  complement  adder  the  "main  adder." 

We  use  several  independent  checkers,  each  of  these  checkers  is  also  an 

ni 

adder,  performing  modulo  2  -I  arithmetic,  where  n^  is  less  than  n.  In 

other  words,  each  checker  is  just  an  n^-bit  long  one’s  complement  arithmetic 

unit.  If  the  n^'s  are  small,  we  assume  that  only  a  small  number  of  errors 

occur  in  the  shorter  adding  units. 

Suppose  we  assume  that  only  single  error  occurs  in  each  short 

adding  unit,  then  instead  of  encoding  the  integers  and  N?  in  the  main 

adder,  we  encode  them  by  different  single-error- correcting  cyclic  AN-codes 

and  add  them  in  the  independent  checkers.  Let  us  assume  that  k  independent 

checkers  are  used,  each  of  length  n^,  l^i£k;  the  single  error  correcting 

ni 

cyclic  AN-code  used  in  the  i-th  checker  is  generated  by  A^,  where  A^B^=2  -1. 

ni 

Since  each  checker  performs  modulo  2  -1  arithmetic,  the  possible 

erroneous  sum  from  the  i-th  checker  is 

ri  =  VRi(Ai(Nl+N2))+V  (5.2.2) 

iijL 

where  R^(x)  denotes  the  residue  of  x  modulo  2  -1.  Since  A^B^=<2  -1, 

(5,2-2)  can  be  rewritten  as 

ri  -  VAi\<W+Eii 


(5.2.3) 


~prrr 


i  <ui*i*  i*" •.*„*',* 
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where  R^  (x)  denotes  the  residue  of  x  modulo  E^. 
i 

Since  generates  a  single  error  correcting  code,  the  single  error 

pattern  can  be  easily  corrected  and  the  residue  (N^+N^)  can  be  obtained 

i 

from  the  output  of  the  decoder.  The  same  argument  holds  for  all  tb*  independent 
checkers,  hence  we  can  determine  all  the  residues  F.  (N.-4-N.)  for  i  =  1,  2,  ...,k 

i.  2 


Suppose  that  the  lengths  of  the  checkers  are  properly  chosen  so  that 
all  the  B^'s  are  pairwise  relctively  prime,  i.e., 


(B^B^)  =  1  for  i^j;  1<  i,  j<  k  (5.2.4) 

where  (a,b)  denotes  the  greatest  common  divisor  (or  g.c.d.)  of  the  integers 
a  and  b.  Writing 


N1  +  N2  “  ^B.  (N1+V  m0d  B1 

N1  +  N2  5  yW  m0dB2 

NL  +  N2  =  Rg  (N^)  mod  Br 


(5.2.5) 


Ly  (5.2.4)  and  the  well-known  Chinese  Remainder  theorem  [49],  we  can  solve  the 

simultaneous  congruence  (5.2.5)  and  uniquely  determine  the  residue  N1  +  N 
k 

modulo  tT  B  .  We  have 
i=i 

k 

N-  +  N_  =  N  +  t  •  rr  B.  ,  for  some  t  (5.2.6) 

1  1  i=l  1 

k 

where  N  =  [(N^+N^)  mod  n  B^  . 

i=l  1 
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Let  us  assume  the  cyclic  AN-code  generated  by  is  of  rate  (rt)^,  and 
the  independent  checkers  are  so  chosen  that  the  product 

2n+2  >  n  B  >  2n+1  (5-2,7) 

i=l  i 


Since  0  N^,  <  2n-l,  we  have 


Nj_  +  N2  <  2(2n-l)  (5.2.8) 


Combining  (5.2.6),  (5.2.7)  and  (5.2.8),  we  have  t  =  0  and 


Nx  +  N2  =  N  (5.2.9) 


Therefore,  the  output  of  the  decoder  gives  the  correct  sum  of  and  N^.  In 

other  words,  instead  of  using  the  longer  n-stage  adder,  th  ;  addition  of  and 

k 

N0  can  be  executeu  by  the  k  independent  shorter  checking  adders;  we  use  E  n. 

1  i=l  1 

elemental  adding  units  to  perform  the  n-bit  binary  addition,  the  redundant 
k 

E  n.  -  n  adding  units  are  incorporated  for  error  correction. 
i-1  1  bi 

Suppose  that  the  B^'s  are  chosen  to  be  B.  a  2  ,  then  the  race  of  the 

code  is 


(rt)£  =  log2  Bi/ni  £  bi/ni  (5.2.10) 


On  the  other  hand,  in  orfcc  to  have  (5.2.7),  the  b^'s  satisfy 

k 

E  b  >.  n  +  1 
i=l  1 


(5.2.11) 


From  (5.2. iO)  and  (5.2.11)  it  is  seen  that 


L  (n  )  5  S  b  /(rt)  £  (n+l)/(rt) 
i_l  i  i  i  max 


- - <  (rt) 

k  \  max 

E  nj 


*>!>  -  (tlogjBj/frOj)) 


(5.2.12 .a) 


n-f2 

k 

£  n. 


£  S  log^l /<«>„*„£  <n+2)/(rt)min 


a  (rt)  . 

min 


(5.2.12.b) 


where  (rt)  =  m®x(rt).  and  (rt)  .  =  m*n(rt).. 

max  i  i  rain  i  i 


The  expression  n/  2  (n.)  is  the  efficiency  of  this  separate  checking  technique. 
i=l  1 

From  (5.2. 12. a)  and  (5.2.12.b),  the  efficiency  is  bounded  by  the  maximum  and 
minimum  rates  of  the  independent  cyclic  AN-codes.  To  increase  the  efficiency 
we  should  choose  the  best  single-error-correcting  AN-code  for  each  independent 
checker.  The  single-error-correcting  An-codes  discovered  by  Perterson  are  all 
perfect,  they  should  be  used  whenever  i t  is  possible.  One  question  is  that  some 
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perfect  codes  are  constructed  with  arithmetic  modulo  2  +1,  while  the  above 

discussion  assumes  that  the  checkers  perform  one's  complement  arithmetic.  How¬ 
ever,  it  is  easy  to  see  that  all  the  above  discussion  also  holds  true  when  some 

ni 

of  the  checkers  perform  arithmetic  modulo  2  +1  ;  this  broadens  the  possibility 

of  using  perfect  codes  for  the  checkers.  Let  us  consider  the  following  example: 

Example  5.2.1:  Suppose  that  a  34-stage  main  adder  is  checked  b'  separate 
checking  technique.  We  use  six  independent  checkers.  The  first  five  checkers 
perform  one’s  complement  arithmetic,  with  lengths  9-bit,  10-bit,  11-bit,  12-bit 
and  15-bit  respectively.  The  cyclic  AN-codes  used  for  chese  five  checkers  are: 


i 

n, 

l 

Ai 

B. 

l 

Dmin 

- 

9 

73 

7 

3 

2 

10 

3.31 

11 

4 

3 

11 

23 

89 

3 

4 

12 

5.7 

3-3  -13 

3 

5 

15 

7.31 

151 

3 

For  the  last  checker,  we  use  a  14-bit  adder  which  performs  arithmetic  modulo 
14 

2  +1.  The  AN-code  for  this  checker  is  a  perfect  sir.gle-error-correcting 

code  (but  not  cyclic)  with  generator  A ,  =  29  and  B,  =  5 '113  =  565. 

o  o 

It  can  be  seen  that  all  the  B^'s  are  pairwise  relatively  prime.  Also, 
6 

TT  B.  =  7 '11  *89 '3  *3  *13 '151 '565 
i=i 

35  34+1 

=  68,  405,  652,  315  >  2  =  2 
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Therefore,  the  34-stage  main  adder  can  be  replaced  by  the  six  independent 


adders. 


The  number  of  redundant  adding  units  is 


£  n  -n  -  (9+10+11+12+14+15) -34  -  37 
i-1  1 

and  the  efficiency  of  the  overall  checked  adding  system  is 


34/71  -  0.48 


which  is  greater  than  one  third. 

The  efficiency  of  the  over-all  checked  system  is  better  than 
the  system  checked  by  the  triplication  method;  also,  not  like  the  tripli¬ 
cated  system  that  all  errors  are  assumed  to  occur  in  one  of  the  three 
identical  processors,  in  separate  checking  approach  we  allow  errors  to 
occur  in  any  independent  checking  unit.  The  assumption  that  only  single 
error  occurs  in  each  independent  checker  is  reasonable  when  the  length 
of  the  checker  is  chosen  to  be  short.  In  addition,  the  decoding  proce¬ 
dure  is  done  separately.  For  each  checker,  the  decoding  is  just  a  single 
error  correction,  which  is  much  simpler  and  much  more  reliable  whan  com¬ 
pared  with  the  approach  of  using  an  over-all  multiple-error-correcting 
AN-code  for  the  main  adder. 


R 
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5.3  Multiple  Iterative  Error-Correcting  Codes  for  High  Speed  Multiplier 
In  a  high  speed  arithmetic  unit,  the  multiplier  is  divided  into 
a  number  of  blocks  with  two  (or  more)  bits  each.  Each  block  is  then  multiplied 
with  the  multiplicand  to  form  partial  sums.  The  partial  sums  are  appropriately 
shifted  and  added  in  a  multi-input  adder  with  minimum  carry  provisions.  If  a 
faulty  circuit  occurs  in  the  multiplier,  trie  error  patterns  have  the  following 
proper ti es: 

a.  Since  the  partial  sums  are  shifted  by  multiples  of  the  length  of  a 
block,  the  erroneous  bits  in  each  partial  sum  will  be  equally  spaced 
when  the  result  is  obtained,  usually  they  span  a  fixed  number  of 
blocks  in  the  binary  representation  of  the  result.  This  error  is 
said  to  be  iterative  in  nature. 

b.  Since  a  stuck-at-1  or  0  type  of  logic  fault  causes  either  a  carry  or 
a  borr'-.rf  error  but  not  both,  all  the  erroneous  bits  of  a  single 
iterative  error  are  of  the  same  polarity.  Utilizing  these  specific 
properties  of  the  iterative  error  patterns,  Chien  and  Hong  [16]  have 
found  a  large  class  of  arithmetic  AN-codes  which  can  correct  any 
iterative  error  caused  by  single  component  failure.  These  codes  are 
good  since  they  possess  high  efficiency  and  can  be  implemented  quite 
easily. 

If  more  than  on  component  of  a  high  speed  multipler  fail,  the  errors 
from  the  output  can  be  recognized  as  several  groups  of  single-iterative-errors. 
We  call  this  type  of  errors  the  multiple-iterative-errors.  Following  the  in¬ 
vestigation  of  Chien  and  Hong's  work,  however,  it  is  difficult  to  generalize 
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their  single-iterative-error-correcting  code.  In  this  section,  we  apply  the 
majority- logic  decodable  AN-codes  constructed  in  the  previous  chapter  to  the 
error-correction  in  high  speed  multipler.  Instead  of  making  correction  at  the 
end  of  multiplication,  we  correct  errors  contained  in  each  partial  sum;  it  can 
be  shown  in  the  following  that  the  majority- logic  decodable  codes  are  suitable 
for  this  purpose,  and  can  correct  any  multiple  iterative-errors. 

First,  let  us  consider  that  the  multiplier  is  divided  into  r  blocks 
and  each  block  contains  m  bits,  then  the  multiplier  is  of  length  rm.  Usually, 
the  multiplicand  requires  the  same  number  of  bits  in  its  binary  representation 
as  the  multiplier,  hence  the  product  requires  at  most  2rm  binary  bits.  A  single 
iterative  error  pattern  E  can  be  expressed  as 

2r-l  . 

E  «  +  2K  E  e,  2m  >  where  0$  k<  m  and  e,  *  0,1  (5.3.1) 

i-0  1 


It  is  easy  to  see  that  a  multiple,  say,  t-iterative-error  Et  can  be 
expressed  as 


E 


t 


2r-l 


£ 

i»0 


'li 


2r-l 

£ 

i=0 


"2i 


2r-l 
£  e 
i“0 


ti 


,mi 


(5.3.2) 

where  0<  k  ,  fc  ,  ...,  k  <  m;  k  's  are  distinct  and  e  *  0,  1  for  all  j’s 

12  t  j  ji 

and  i'a. 

The  error  patterns  (5.3.1)  or  (5.3.2)  are  presented  in  the  final 
result  at  the  output  of  a  high  speed  multiplier.  If  we  consider  the  partial 
error,  which  is  the  error  contained  in  each  partial  sum,  we  have  the  following 


form 
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partial  error  at  the  i-th  partial  sum 


k  +mi 
2  1 


+  e 


2i 


k  -hni 
2  L  + 


+  e 


ti 


k  -Hni 
2  t 


(5.3.3) 


where  e^  *  0,  1  or  -1.  0^  k^,  ...  <  m  and  all  the  k^*s  are  distinct.  There¬ 

fore,  each  partial  error  of  a  t-iterative  error  pattern  can  be  considered  as  a 
random  error  of  arithmetic  weight  t  or  less.  This  suggests  a  possible  way  of 
correcting  any  t-iterative-errors . 

We  encode  the  multiplicand  by  a  t- random-error-correcting  AN-code. 
Instead  of  making  correction  of  the  final  product  we  use  the  t-error-correcting 
code  to  correct  the  partial  errors  contained  in  the  partial  sum.  Certainly,  the 
decoding  of  the  t-error  code  must  be  very  easy,  otherwise  the  approach  is  highly 
impractical.  The  majority- logic  decodable  AN-codes  constructed  in  Chapters  3  and 
4  immediately  render  the  potential  of  being  the  t-random-error-correcting  codes. 

Design  of  the  Codes 

Suppose  that  both  the  multiplier  and  the  multiplicand  require  rm  bits 

in  their  binary  representations,  then  the  product  needs  2rm  binary  bits.  We 

want  to  design  an  arithmetic  AN-code  of  length  2rm,  with  the  generator  A  a  pro- 
2rtn 

per  divisor  of  2  -1,  then  the  code  is  cyclic  and  all  the  code  words  AN  form  a 

2rtn 

principal  ideal  in  the  ring  of  integers  modulo  2  -1. 

Let  8  be  a  divisor  of  rm,  then  rm  »  d*s,  for  some  d.  The  integer 
2^  l)8  +  2^  +  ....  +  2®  +  1  can  be  chosen  as  the  6enerator  of  the  code. 
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First,  we  see  that 

A  -  2W-Ds  +  2(d-2)s  +  +  2s  +  1 

-  (2rm-l)/(2s-l)  (5.3.4) 

Since  aJ (2  ~1),  the  integer  2  -1  is  a  code  word  of  weight  2  and  the  code 
generated  by  A  is  of  minimum  distance  only  2.  It  seems  that  the  code  has  no 
error  correcting  capability,  yet  it  can  correct  all  t-iterative-errors  if 
2t  +  IS  d. 

Let  the  multiplicand  be  encoded  as  a  code  word  AM,  and  let  the 
multiplier  be  N.  In  a  high  speed  multiplier  the  binary  form  of  N  is  divided 
into  r  blocks,  each  block  contains  m  bits;  hence  N  can  be  expressed  as 

r-1 

mi 

N  *  £  n  2  ,  where  0<  n  <  2m-l  (5.3.5) 

i=0  1  1 

Denoting  as  the  i-th  partial  sum,  we  have 
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where  Y1  is  the  erroneous  result  at  the  i-th  stage  of  multiplication.  Denoting 


the  residue  of  x  modulo  2m-l  by  R.m(x),  (5.3.8)  can  be  written  as 


Rrm<V  "  «tJ*f<^i.i)  +  Epil 


(5.3.9) 


If  t  or  less  iterative  errors  occur  in  the  multiplication,  then  from  (5.3.3) 


the  weight  of  Epi  is  less  than  or  equal  to  t.  Since  A  *  (2rm-l)/(2S-l)  =  2^d_1^s 


+  ...  +2  +  1  and  ds  =  rm,  the  residue  of  the  integer  AMnt  x  modulo  2rm-l, 

Rrm^**ni-1^  contalns  d  indentical  blocks  in  its  binary  representation  form, 

with  each  block  df  length  s.  From  (5*3.9),  it  is  seen  that  OS  |e  I  <  2rm-l 

pi 1  * 

W(Epi)  <  t;  hence  the  majority- logic  decoding  scheme  developed  in  Chapter  4  is 

applicable  to  the  residue  R^O^).  A  majority  decision  on  the  binary  bits  of 

Rrm^Yi^  can  correctly  yield  the  partial  error  pattern  Epi  provided  that 

d  2  2t  +  1.  The  same  majority- logic  decoding  scheme  can  be  applied  to  all 

Rrm^?  for  1  *  1»  2,  ...,  r,  and  all  the  partial  errors  E  , ,  E  _ ,  ...  E 

pl  p2  pr 

can  be  determined. 

The  design  of  the  above  code  is  based  on  the  assumption  that  both 
the  multiplicand  and  the  multiplier  have  rm  bits  in  their  binary  forms.  As 
a  matter  of  fact,  this  result  can  be  easily  generalized  to  any  case  when  the 
multiplicand  and  the  multiplier  are  of  different  length,  say,  r^  bits  for  the 
multiplicand  and  r2m  bits  for  the  multiplier.  In  these  cases,  the  generator 
of  the  code  should  be  modified  as 

*'  -  +  2M-2>.  +  ....  +  2»  +  ! 

where  ds  *  rjm,  and  the  length  of  the  code  is  (r.+r  )ra.  The  generator  A'  may 

(ri+Om  1  L 

not  be  a  divisor  of  2  -1  and  the  code  may  not  be  cyclic.  Nevertheless, 


9' 


the  error  correcting  procedure  of  the  code  is  still  the  same. 

Also,  it  should  be  remarked  that  the  code  devised  above  can  correct 
iterative  errors  caused  by  any  possible  logic  faults  (not  only  the  stuck-at 
type  of  fault)  in  the  circuitry  of  a  high  speed  multiplier.  The  erroneous  bits 
of  a  correctable  iterative  error  pattern  are  not  restricted  to  have  the  same 
polarity. 


Rate  of  Code 

The  information  rate  (or  the  efficiency)  of  an  arithmetic  code  is 

defined  as 


Rate  *  1  -  [  (.lumber  of  bits  in  the  binary  form  of  A)/length  of  code] 

(5.3.10) 

For  the  code  of  length  2rm,  generated  by  A  *  (2rm-l)/{2*-l) ,  where  ds  *  rm, 
the  rate  can  be  readily  calculated  as 


Rate  (t-iterative-error-correcting  code) 


rm  +  (rm-ds+s) 
2rm 


(5.3.11) 


It  is  interesting  to  see  that  the  rate  of  any  t-iterative  error 
correcting  code  is  greater  than  0.5.  For  t  »  1  (single-iterative-error),  we 
take  d  ■«  2t  +  1  *  3.  Then  3s  *  rm  and  the  rate  is 
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Rate  (single-iterative-error-correcting  code' 

1  .  i  1.1.1  2 

2  2m  2  2  3  3 

Therefore,  the  rate  of  any  t-iterative-error-correcting  code  of  length  2rm  lies 
in  the  range 

■r  +  —  S  rate  (t-iterative-error  code)  S  \  (5.3.12) 

z  m  j 

For  different-length  multiplicand  and  multiplier,  there  may  exist  some 
t-iterative-error-correcting  codes  with  efficiency  even  better  than  2/3. 

From  the  above  discussion,  we  see  that  the  construction  of  the  code 
is  simple,  the  efficiency  of  the  code  is  satisfactory  and  more  importantly,  the 
decoding  can  be  achieved  by  using  simple  one- level  majority  logic.  On  the  other 
hand,  since  the  error  correcting  is  done  at  each  stage  of  the  multiplication 
when  a  partial  sum  is  obtained,  there  is  definitely  a  time  delay  in  the  soeed 
of  multiplication. 
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6.  CONCLUSIONS 


6.1  Conclusions 

The  goals  of  this  thesis  have  been  to  investigate  the  two  unsolved 
problems  in  the  theory  of  arithmetic  coding  mentioned  in  the  introductory 
chapter,  to  develop  new  checking  techniques  for  computer  arithmetic  and  to 
examine  the  theoretical  structure  of  the  arithmetic  codes.  For  the  code  con¬ 
struction  problem,  several  new  classes  of  multiple-error-correcting  arithmetic 
codes  are  found.  By  using  some  interesting  properties  of  the  binary  representa¬ 
tion  forms  of  certain  integers  the  error  control  capability  of  the  codes  can 
be  determined  analytically.  The  synthesis  of  these  classes  of  codes  are 

simple,  in  fact,  for  any  two  relatively  prime  integers  n^  and  n^,  the  generator 

tun 

of  the  code  is  chosen  to  be  (2  -l)/(2  -1)(2  -1).  The  fact  that  pj^-code, 

which  is  a  subclass  of  the  codes  constructed,  possesses  a  root-distance  relation¬ 
ship  similar  to  that  of  the  BCH  bound  indicates  additional  analogy  between  arith¬ 
metic  AN-codes  and  cyclic  parity-check  codes. 

For  the  decoding  implementation  problem,  the  majority-logic  decoding 
scheme  for  arithmetic  codes  is  devised.  This  is  a  completely  different  approach 
from  the  attempts  made  by  using  various  search  techniques  [29, 30] ,  the  impleb 
mentation  of  the  codes  can  be  achieved  by  using  majority  logic  alone.  It  has 
been  shown  that  the  arithmetic  AN-codes  constructed  systematically  in  this  thesis 
are  majority- logic  decodable.  The  majority-logic  decoding  scheme  is  applied  to 
other  classes  of  AN-codes  and  a  generalized  multi-step  decoding  algorithm  is 
obtained.  The  number  of  steps  of  majority-logic  decision  required  in  a  multi- 
step  decoding  can  be  shown  to  equal  the  number  of  distinct  prime  factors 
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conatained  in  the  block  length  of  the  code.  It  has  also  been  noticed  that  the 
minimum  distance  of  these  majority- logic  decodable  codes  can  all  be  determined 
exactly. 

A  new  checking  technique  for  the  binary  adders  has  been  developed. 

This  separate  checking  technique  provides  a  more  reliable  arithmetic  error 
checking  model,  the  redundant  circuits  used  are  less  than  that  required  for  a 
triplicated  system.  Since  the  arichmc.tl c  coding  scheme  is  used  separately  for 
each  shorter  adding  unit,  the  single- error  assumption  is  made  and  the  decoding 
is  much  simpler.  Furthermore,  the  application  of  majority- logic  decodable 
arithmetic  AN-codes  to  the  errot  -orrection  in  high  speed  multiplier  has  been 
examined. 

Some  theoretical  results  about  the  structure  of  arithmetic  AN-codes 
have  been  obtained.  By  using  the  concepts  of  number  theory,  the  upper  and  lower 
bounds  on  any  cyclic  AN-code  have  been  formulated;  also,  a  fast  algorithm  for 
finding  the  binary  representation  forms  of  certain  integers  has  been  derived. 

In  summary,  the  results  presented  in  this  thesis  gives  a  better  belief  that  the 
use  of  arithmetic  coding  for  the  error  control  in  digital  computing  systems  is 
indeed  very  promising. 

6.2  Suggestion  for  Future  W Prl'. 

It  seems  that  the  next  major  advance  in  arithmetic  coding  from  the 
theoretical  viewpoint  will  be  the  complete  formulation  of  the  class  of  cyclic 
AN-codes  analgous  to  the  Bose-Chaudhuri-Hocquenghera  parity-check  codes.  When 
this  problem  is  solved,  the  computer  designer  interested  in  incorporating 


..  ..x*. 
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coding  for  arithmetic  error  control  will  have  an  easy  way  of  constructing  codes 
From  practical  point  of  view,  the  separate  checking  technique  discussed  in  this 
thesis  will  find  its  way  into  actual  computer  hardware  designed  to  operate 
reliably  either  faster  or  for  a  longer  time  than  non-coded  equipment.  How  to 
syste-nactically  choose  the  various  length  of  the  independent  checking  unit 
so  that  the  redundancy  is  minimized  would  be  a  problem  which  requires  further 
examination. 

The  implementation  of  decoder  by  using  majority- logic  increases  the 
possibility  of  widespread  application  of  arithmetic  coding.  Future  work  will 
be  the  discovery  of  other  more  efficient  majority-logic  decodable  arithmetic 
AN-codes  and  other  practical  decoding  methods  based  on  the  concepts  of 
majority-logic  decision. 
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